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We extend earlier results on the relation between the dimensionless tunneling channel conduc- 
tance g and the fractional Coulomb blockade peak splitting / for two electrostatically equivalent 
dots connected by an arbitrary number iV c h of tunneling channels with bandwidths W much larger 
than the two-dot differential charging energy Ui. By calculating / through second order in g in 
the limit of weak coupling (g — > 0), we illuminate the difference in behavior of the large-iVch and 
small- A^h regimes and make more plausible extrapolation to the strong-coupling (g — » 1) limit. For 
the special case of N c h = 2 and strong coupling, we eliminate an apparent ultraviolet divergence 
and obtain the next leading term of an expansion in (1 — g). We show that the results we calculate 
are independent of such band structure details as the fraction of occupied fermionic single-particle 
states in the weak-coupling theory and the nature of the cut-off in the bosonized strong-coupling 
theory. The results agree with calculations for metallic junctions in the AT c h — > oo limit and improve 
the previous good agreement with recent two-channel experiments. 

PACS: 73.20.Dx,71.45.-d,73.40.Gk 
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I. INTRODUCTION 



The opening of tunneling channels between two quan- 
tum dots leads to-, an erosion of the individual dots' 
Coulomb blockade.^ For a pair of electrostatically identi- 
cal quantum dots (see Fig. 1 for a schematic view of the 
double-dot structure), the progress of this erosion can 
be chronicled by tracking the splitting of the Coulomb 
blockade conductance peaks as they evolve from doubly 
degenerate single-dot conductance resonances to that of 
nondegenecate double-dot peaks with twice the original 
periodicityJaQ For a system in which the tunneling chan- 
nels can be treated as having the same individual conduc- 
tances and in which the Coulomb charging energies are 
large compared to the single-particle level spacings but 
small compared to the tunneling channel bandwidths, the 
fractional peak splitting / can be expressed as a func- 
tion of two parameters: A7 ch , the number of tunneling 
channels between the two dots, and g, the dimensionless 
conductance per tunneling channel. (In this paper, the 
conductances indicated are always dimensionless conduc- 
tances, by which we mean the actual conductance divided 
by the conductance quantum, e 2 /h.) 

In particular, for weakly coupled dots (g — > 0), the 
fractional peak splitting can be expressed perturbatively 
as a sum of terms of the form a mtn {N c \ l ) m g n , where 
1 < m < n i-P-ffl a m ,n is independent of N c h and g. 
Previous workuTQ has produced the leading term in this 
expansion. However, as this term is simply linear in 
the total interdot tunneling conductance, gtot = N c ^g, 
it does not effectively distinguish between behavior in 
the large-A/^h and small- AT c h limits. To make such a 
distinction, one must calculate to second order in g, in 
which case one obtains two sets of terms, one set propor- 
tional to N c hg 2 = gt ot /Nch and the other proportional to 



(N c h) 2 g 2 = 3tot- The (A^h) 2 terms should agree with re- 
sults from the theory of metallic junctions, in which the. 
leading terms in a large- N c h expansion are calculatedE-El 
In the strong-coupling limit (g —* 1), a dramatic de- 
pendence of the peak splitting on A^h has already been 
found. In the cases of N c h = 1 and AT c h — 2, the leading 
nontrivial terms have been found to be or/anortional to 
y/1 — g and (1 — ffUfttt — g), respectivelyjatffl and it has 
been hypothesizeofcltLil that, for N c h > 2 but finite, the 
leading nontrivial term is proportional to (1 — g) Nch ^ 2 . 
This last suggestion appears—to correspond to calcula- 
tions in the large- N c ^ limit ,113113 where the effective charg- 
ing energy has been found to be proportional to e~ fftot / 2 , 
which is equivalent to (1 — gtot/^Vch)^ 011 ^ 2 in the limit 
N ch — > oo. 

Despite the recent progress in study of the strong- 
coupling lirait r _for the case of most direct experimen- 
tal interest J3ut2l N c h — 2, the leading-term calculation 
fails to be completely satisfactory for at least three rea- 
sons. The first is that this calculation does not answer the 
question of whether the coefficient of (1 — g) ln(l — g) 
is affected by the manner in which the ultraviolet cui« 
off is imposed in the low-energy Luttinger liquid theory.Q 
The second is that the coefficient of the sub-leading terra 
linear in (1 — g) is both unknown and naively infinite.Ll 
Finally, there is the worry — which also applies to the 
weak-coupling result — that, for A^ c h = 2, interpolation 
between the solutions for weak- and strong-coupling is 
difficult because the respective /-versus-g curves do not 
come especially close.Q 

This paper addresses these three concerns for the two- 
channel problem and also extends earlier results for the 
general A^h-channel problem in the limit of weak cou- 
pling. In so doing, it illuminates the difference between 
large-A^h and small-A^h behavior for g 0, creates a 
theory that can be more realistically compared to exper- 



1 



imental results for AT c h = 2, and argues for the univer- 
sality of the results, which should be independent of the 
nature and magnitude of the ultraviolet cut-offs. Sec- 
tion II presents the g 2 extension of the weak-coupling 
theory and checks the result against calculations in the 
-/V c h — > oo limit. Section III gives the (1 — g) correction 
to the leading dependence in the strong-coupling limit 
for iV c h = 2 and shows a plot of the experimental results 
and revised theoretical predictions for two-channel inter- 
dot junctions. Section IV argues that the strong coupling 
results of Section III are independent of the nature of the 
way the ultraviolet cut-off is imposed and do not change 
when one allows the fermionic theory to stray from half- 
filling. Section V summarizes the results, and Appen- 
dices A and B present technical details of calculations in 
Sections II and III, respectively. 




FIG. 1. Schematic diagram for the double dot. Negative 
potentials are applied to each of the gates to form the dou- 
ble-dot structure. The gate potentials V g i and V g 2 control 
the average numbers of electrons on the dots. These are 
the potentials that are varied to see the Coulomb block- 
ade. Vb controls the rate of tunneling between the dots. 
Vxi and Vx2 control the rate of tunneling to the adjacent 
bulk two-dimensional electron-gas (2DEG) leads. For cal- 
culations of the double-dot energy shifts, tunneling to the 
leads is assumed negligible compared to tunneling between 
the two dots. In measuring the total channel conductance 
Gtot , however, the potentials Vxi are turned off so that each 
dot is strongly connected to its lead. The side-wall poten- 
tials V s i and V S 2 are fixed. 



II. THE WEAK-COUPLING LIMIT FOR 
ARBITRARY N C n 



For weakly coupled quantum dots, we use a model 
"site-to-site" hopping Hamiltoniano and calculate per- 
turbatively in the tunneling term Ht- 



Hr — ^ ^ ^ ek^iko- , 
i=l a k 

H c = U 2 {h-p/2f, 



(1) 



a kik 2 



As in Ref. 4, in these equations, i is the dot index; a 
is the channel index; k is the index for all internal de- 
grees of freedom not included in the channel index; Hq 
is the part of the electrostatic potential energy that is 
affected by interdot tunneling; n is half the difference in 
dot occupation numbers, n = (fi2 — ni)/2; p is a differ- 
ential gate voltage parameter and is restricted to values 
between and I (as permitted by the system's unit peri- 
odicity); and U 2 is the differential charging energy, which, 
for electrostatically equivalent dots, is given by the for- 
mula U2 = e 2 /(Cs + 2Ci n t), where C[ nt is the interdot 
capacitance and Cs is the total single-dot capacitance 
minus the interdot capacitance. If the dots are not elec- 
trostatically equivalent, the formula for U 2 and the defi- 
nition of p are more complicated. uu However, the model 
is still applicable, and the results for f p can still be used 
to obtain the peak splitting. 

These calculations are made palatable by assuming 
that U2 is much smaller than the tunneling-channcl band- 
width W yet much greater than the average intrachanncl 
level-spacing 5: W 3> U2 S> S. This assumption leaves 
us with a theory that we can consider to be in the con- 
tinuum limit and that we can hope to be independent 
of ultraviolet cut-offs. As the bandwidth is presumably 
of the order of the Fermi energy ep, these assumptions 
are reasonable for the micrometer-sized dots of Waugh 
et aZ.,Bfl for which e F ~ 10 meV, U 2 « 400 ^eV, and 
S w 30 fieV. 

As in Ref. 4, our primary goal is to calculate the frac- 
tional peak splitting / — i.e., the ratio of the distance be- 
tween split Coulomb blockade subpeaks for a given g and 
their maximal separation in the strong-coupling (g — > 1) 
limit. It was shown in Ref. 4 that, if the total number of 
electrons on the two dots is assumed even, the problem 
of solving for / is a corollary to the problem of solving 
for a more general quantity f p , which characterizes the 
ground-state energy of the double-dot when the differ- 
ence between the external potentials applied to the dots 
is nontrivial and the total number of electrons on the two 
dots is fixed and even. Recall the equation for f p : 



A - A p 

U 2 p 2 /4: 



(2) 



where A p is the shift in the ground state energy induced 
by tunneling at a given value of the gate voltage param- 
eter p and U2P 2 /4 is the difference between the unper- 
turbed ground state energies for the given p and p = 0. 
In Ref. 4, it was shown that, for symmetric dots, 



H =H K + H C 



f — fp=i 



(3) 
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In the same work, it was determined that f p exhibits the 
following leading behavior as g — > 0: 



4 X) = K 1 - P) Ml ~ p) + (1 + p) ln(l + p) 



TT^p^ 



(4) 



where ^ = W/U2 3> 1. Thus, the corresponding leading 
behavior for / is 



/W^jv c , 9 + o( 9 M s 2 ) 



(5) 



Extending perturbation theory beyond this result — 
i.e., beyond first-order in g — requires some laborious 
computation. The next-leading contributions come from 

two sources. The first, which we shall call fj> 2A \ arises 
from a combination of the second-order energy shift that 
has already been calculated and the second term in the 
formula that relates the tunneling amplitude t to the 
channel conductance g. [The first term in this formula 
was used to derive Eq. (Q).] The second source of g 2 

terms, fp 2B \ is the shift in the ground-state energy pro- 
vided by terms that are fourth-order in t. 

The first contribution is relatively easy to calculate. 
The equation for g in terms of t has been derived for 
half-filling in Ref. 14 and for arbitrary filling in Ref. 4. 
In the latter calculation, the system is assumed to have 
a constant density of states between single- particle en- 
ergies e and (e -I- W), the density of states being zero 
elsewhere. The system's level of "filling" is then charac- 
terized by the filling fraction F = (ep — cq) /W, where cf 
is the Fermi energy. In accordance with the half-filling 
result, one then finds the following: 



9 = 



4x 



|1 + (1 + ir)) 2 x\' 



(6) 



where x = (^t/S) 2 and 77 = (I/71") ln[F/(l— F)]. Inverting 
this expression, one discovers that 



^2 



9 
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T 



9 + 0(g 2 ) 



(7) 



Consequently our first g 2 term is equal to the right side 
of Equation (0) multiplied by (1 — i] 2 )g/2: 



= (l-n 2 ) Nch92 



[(1 -p) ln(l-p) 



2ttV 

+ (l + p)ln(l + p) + <rVAA)]. (8) 



This term is of the expected form a^ 2 ^ N^g 2 , where 
is a function of p. 

On the other hand, a^ 2 ^ is dependent on the filling 
fraction F, a fact which appears to imperil our dreams of 
a theory that is universal in that it is insensitive to the 
details of the high-energy behavior (including whether, 
for example, certain high-energy states exist and there- 
fore have a role in determining the filling fraction F). 



,(2A) 



We shall see, however, that the F-dependence of fp 2 ^ 
actually serves our end, for it exactly cancels the in- 
dependence of fp 2B ^ . As a result, we can further conclude 
that, through second order in the channel conductance g, 
expression of the fractional peak splitting in terms of the 
channel conductance is not only convenient for compari- 
son with experiment but is also necessary and sufficient 
for constructing a result that can be hoped to be univer- 
sal. 

To support this claim, we must actually determine the 
value of f p . Sadly, it cannot be obtained as effortlessly 

(2A) 

as fp . There are 24 separate terms that contribute 
to the fourth-order energy shift. One 12-member subset 
consists of terms proportional to (iV c h) 2 ; the second con- 
sists of those simply linear in N^. All but four of the 
twenty-four terms correspond to a specific series of four 
tunneling events that begin and end with the double- 
dot system's unperturbed ground state. The remaining 
four, which belong to the (iV c h) 2 subset, correspond to 
the fourth-order terms in Rayleigh-Schrodinger perturba- 
tion theory that are products of the second-order energy 
shift and a propagator squared. These four have, been 
described by Grabert as diagrams with insertions.u 

In general, the nature of the twenty-four fourth-order 
terms is most digestibly summarized via a diagrammatic 
representation that looks essentially like one of time- 
ordered single-particle diagrams (see Fig. 2). Despite the 
superficial single-particle nature of this representation, it 
is important to remember that the propagators that enter 
into the energy calculations are the propagators for the 
entire double-dot system, which depend upon both the 
tunneling particles' individual kinetic energies and the 
system's multiparticle potential energy. The presence of 
the multiparticle potential energy makes it impossible to 
reduce the calculation to the normal Feynman diagrams, 
for which one can write the problem entirely in terms 
of single-particle propagators. The presence of exchange 
terms, which do not appear among the diagrams propor- 
tional to (iV c h) 2 , makes a pseudo-single-particle represen- 
tation necessary. 

Within this time-ordered perturbation theory scheme, 
the individual fourth-order terms are plagued by both 
ultraviolet and infrared divergences. Every term is di- 
vergent as the bandwidth goes to infinity and four of the 
(A^ch) 2 terms are divergent as p — ► 1. (A different set of 
four is divergent as p — > — 1.) From the result for fp X \ 
we might hope to cancel the ultraviolet divergences and 
to obtain an answer for the ground-state energy that is 
infrared-singular but not infrared-divergent. Indeed, as 
Grabert has notedfj the ultraviolet divergences of the 
(A c h) 2 terms must drop out since, in the limit U2 — > 0, 
these terms correspond to disconnected diagrams or in- 
sertion diagrams that exactly cancel one another and 
thus do not appear as Feynman diagrams. In contrast, 
the A c h diagrams do have nontrivial Feynman-diagram 
analogs. As a whole, they correspond to a single totemic 
Feynman diagram — an individual ring marked by four 
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FIG. 2. Diagrams for half of the (a) fourth-order, (N^) 2 
terms and (b) fourth-order, N c h terms. The remaining 
terms are represented by diagrams that are mirror images 
of these. A vertical dashed line is drawn for each of the 
m particles that tunnels at least once from one dot to the 
other. This line stands for the corresponding particle's ini- 
tial state, a state that must be filled at the end of the 
four tunneling events in order to recover the unperturbed 
ground state from which the system starts. A particle be- 
gins at the bottom of its vertical initial-state line. Particles 
in dot 1 propagate upward and rightward. Particles in dot 
2 propagate upward and leftward. A tunneling event for 
a particle is signaled by a solid dot that coincides with a 
bend in the particle-propagation path. Each particle must 
end on one of the dashed vertical lines, meaning that it 
ends in the single-particle state that corresponds to that 
line. Insertions (see Sec. II) are represented by triangles 
that project off a single-particle propagation line. If the 
projection points up, the insertion corresponds to the term 
in the second-order energy shift for which a particle tun- 
nels off the dot occupied by the propagating particle. If 
the projection points down, the insertion corresponds to 
the second-order term for which a particle tunnels onto the 
dot occupied by the propagating particle. In the absence of 
exchange, all particles end on their own initial-state lines. 
A two-particle exchange carries a minus sign and results in 
each of two particles ending on the other's initial-state line. 
Three-particle exchange carries no sign (alternatively, one 
can view it as carrying two canceling minus signs) and re- 
sults in each of three particles ending on one of the others' 
initial-state lines. 



tunneling events. The ultraviolet divergences of these 
diagrams are therefore expected to be persistent but ir- 
relevant because we are interested only in the relative 
shift between the ground-state energies for arbitrary p 
and for p = (recall Eq. ||). Accordingly, we expect 
that, when one subtracts the fourth-order shift for p = 
from that for arbitrary p, the fourth-order terms produce 
a result that is neither ultraviolet- nor infrared-divergent 
but is infrared-singular as \p\ — > 1. A brief summary of 
the actual calculation of these terms follows. Those in- 
terested in more detail are invited to peruse Appendix A, 
which offers a fuller description of the calculation of the 
(-/Vch) 2 diagrams and a step-by-step computation of the 
contribution from one representative iV c h term. 

For the less scrupulous, there are still a few facts wor- 
thy of note. A prominent feature of the fourth-order cal- 
culation is that each term involves the integration over 
four energy variables (ej, where i ranges from 1 to 4) of 
a product of three propagators. In the (JV c h) 2 diagrams, 
the energy variables "pair off" : e\ and £3 only appear as 
parts of the combination e\ = {e% — ei), and £2 and €4 
only appear as parts of the combination en = (64 — 62). 
As a result, calculation of these terms reduces to the per- 
formance of double integrations over t\ and en — albeit 
with a nontrivial density of states. 

The N c h diagrams cannot be handled in this way, for 
they involve particle exchanges that frustrate any desire 
to pair off the energy variables. The quadruple integra- 
tion over the cannot be eluded. It can, however, be 
expedited by differentiating twice with respect to p while 
integrating out the energy variables and, then, integrat- 
ing twice with respect to p in the end. One might worry 
about the fact that, by differentiating twice with respect 
to p, one has lost knowledge of terms constant and linear 
in p. However, these terms are unimportant. As noted 
in Ref. 4, the ground-state energy (perturbed or unper- 
turbed) is symmetric in p. Therefore, terms linear in p 
must cancel out of the fourth-order energy shift when 
all the terms are summed. Constant terms are similarly 
negligible since, as usual, we are only interested in the 
relative energy shift A p — Ao. 

After the aforementioned tricks for calculating the 
(^Vch) 2 and A^h diagrams have been used, the only real 
wrinkles that remain are integrals of the form 

Jo {x + A) 

where < \A\ < B, R is either F or (1 — F), and, as 
before, ip = WjUi. The symbol V indicates that, for 
A < 0, only the principal value of the integral is cal- 
culated. These integrals can be done by rewriting the 
argument of In (2; + B) as [(a; + A) + (B — A)] and Taylor- 
expanding about (x + A) for (B — A) < \x + A\ and about 
(B — A) for (B — A) > \x + A\. The result of such an inte- 
gration may be sensitive to whether the system is below 
half-filling [F < (1 - F)], at half-filling [F = (1 - F)}, 
or above half-filling [F > (1 — F)]. However, the system 
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as a whole has particle-hole symmetry, so one expects 
that the final result — once all the terms are summed — is 
symmetric under exchange of F and (1 — F), If there is 
no jump discontinuity when the system is precisely half 
full, the result for F < (1 — F) should determine the 
answer for all "finite" F, by which we mean all F such 
that Ftp, (1 — F)tp ^> 1. This thesis has been explicitly 
confirmed. 

Indeed, the (iV c h) 2 part of the fourth-order relative en- 
ergy shift is found to be independent of the filling frac- 
tion. Its contribution to f p has a rather lengthy explicit 
form: 



{-— p + 4(1 - p) ln(l - p) 



P,(Wch) 2 



4ttV 6 



I-P 2 



ln 2 (l-p)-2(2-p)ln[2(2-p)] 



+ i±£-ln(l + p)ln(l-p) 

- 2 ln(3 - p) ln(l - p) + ln 2 (3 - p) 

+ (3 - p) 2 (1 - p) Ml- P )-ln(3-p)] a 
(_l)n+l /l_p> 



(5-4p + p 2 )£ 



3-p 



5 In 2 3 



8 In 2 + 5k 



+ [p -> -p] } , 



(9) 



where the contents of the last pair of brackets indicate 
that one sums over all the terms in the curly braces again 
after replacing p with —p and the quantity k is given by 



oo 

£ 

n=l 



(-1) 



n+1 



(10) 



It has been confirmed numerically that the derivative 
of /p(jv h )2 agrees with a calculation by Grabert in the 
N c h — > oo limit Jl More precisely, the derivative of this 
term with respect to p corresponds exactly to Grabert 's 
calculation of the leading corrections to the ground-state 
expectation value of h. The value at p — of the second 
derivative of this term with respect to p (and of the first 
derivative with respect to p of Grabert 's n correction) 
should agree with the large- N c b "effective capacitance" 
of Golubev and Zaikin.Q However, this is found not to be 
so, the g 2 correction calculated by Golubev and Zaikin 
being nearly a factor of 20 smaller than that computed 
here. The present calculation therefore provides an im- 
portant check on the large- N c \ x calculations, resolving an 
apparent contradiction in the literature. 

There are no corroborating calculations for the terms 
that are linear in N c ^ as these are negligible in the large- 
-/V c h limit. However, knowing that / = 1 when g = 1 
and that (0.14)A^ c hg, one might conjecture that 

the sign of the g 2 term changes from positive to negative 



when iV c h is of order 10. With respect to the expansion 
of /, such a crossover would imply that the coefficient of 
Nchg 2 is positive and approximately 10 times the size of 
the negative coefficient of (iV c h) 2 g 2 . 

To check this conjecture, we need to know the value of 
the fourth-order, linear-in-iV c h contribution to the frac- 
tional peak splitting. Our results for this quantity are 

fOSL = S^VCl - p) ln(l - P) + |(1 - P) ln 3 (l - p) 
-2(2 + ln2)(l-p) ln 2 (l-p) 

+ 4 ^2 + In 2 - In 2 In 3 - k + y ^ (1 - p) ln(l - p) 
-2 ^y - 2k- 21n21n3^ [(3 - p) ln(3 - p) -31n3] 

- H[( 3 - p) In 3 (3 - p) - 3(3 - p) In 2 (3 - p) 

+ 6(3 - p) ln(3 - p) - 3 In 3 3 + 9 In 2 3-18 In 3] 

- 2 In 2[(3 - p) In 2 (3 - p) - 2(3 - p) ln(3 - p) 

-31n 2 3 + 61n3] 
n 

-2£>(P) 

i=l 

+ [p -"/>]}, (11) 



where k is given by Eq. (|10|), r/ = (1/tt) ln[F/(l - F)], 
and the Ai(p) are defined below: 



Mp) 



P d ^° _ X ^ ^ ^ ~~ X ^ 
X l-x 



M p)=2jyj p ~ x ^i-_f^-^ 

A 3( p)=-21n 2 r <fa (P-»)M3-») 
Jo 1 ~ a: 

A 4 (p) - 2(1 - p) ^"dx In In 

1 1 

x 



3 — x 1 — x / 

A 5 (p)=2 /^(p-^lnfl^) (^--.J- 
Jo \3 — x/\3 — a; 1 — a; 

^.-J'd^-'"'!-'""-' 

Jo 3 x 

i / \ „, „ f , (p - ar) ln(l - i) 

A 8 (/o) = -21n2 / dx^- ^— i i 

Jo s — x 

rP 



A 9 (p) = -2(3- p) f dx In 
Jo 



2 — x\ , ( 3 — p 
In 



3 — x 
I 1 



3 — x 



j 4io(/°) = — 2 / dx(p — x)hx 



3 — a; 1 — a; 

2 - x 



3 — a; 
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*(—-—) 

\3 — x I — x J 

rP 

An{p) = - l dx ln 2 (l -x)ln(l +x) (12) 
Jo 

The characterization of the fourth-order energy shift 
is now essentially complete. The result is more unwieldy 
than we would like. However, there are a few highlights 
that are easy to draw out. As expected, the fourth-order 
shift is neither ultraviolet- nor infrared-divergent but is 
singular as \p\ — > 1, the leading singularities being in 
agreemertt.with an earlier calculation by Glazman and 
Matveev.Ej In addition and quite gratifyingly, the solu- 
tion is independent of the filling fraction F. As discussed 
earlier, the dependence of fp 2B ^ on the filling fraction, 
which is concentrated in the rj 2 term of the first line of 
Eq. (pl|), exactly cancels that of Eq. (g). Hence, there 
is some reason to believe that, when expressed in terms 
of the channel conductance g, the result is universal in 
the sense that it is independent of the details of the band 
structure for energies much greater than U2, where U2 is 
much less than the bandwidth W. 

It is difficult to get a better handle on this algebraic 
smorgasbord by mere inspection. One can add some pre- 
cision to the picture of what has been accomplished by 
first assembling the g 2 terms of f p and then plugging in 
p = 1 to obtain the contribution to the symmetric-dot 
fractional peak splitting /. Upon recalling that 

f (2) _ f(2A) , f(2B) A2B) 

one can evaluate the A\ integrals numerically for p = 1 
to obtain 

/< 2 > a [0.1491] N ch g 2 - [0.009798] {N ch ) 2 g 2 . (13) 

We see that the conjecture about the (A c h) 2 and A c h 
contributions to is correct: the terms have opposite 
sign, and the ratio of their magnitudes is on the order of 
10. For the case of A c h = 2, the g 2 term provides the 
desired upward correction to the /-versus-g curve. 

Before specializing to the result for A c h = 2, we should 
explore the consequences of having a term proportional to 
Nchg 2 . This term makes the result sensitive to the "fine 
structure" of the interdot conductance. As remarked 
in the introduction, terms of the form (N^g) 71 can be 
rewritten as a simple power of the total conductance be- 
tween the dots: {N^g) n — (<?tot)"- Should the conduc- 
tances in the various tunneling channels be allowed to 
differ, the form of these terms when written in terms of 
5tot would remain unchanged. The only alteration would 
be in the equation for <7 to t itself, which would revert to 
the more fundamental form 

fftot = X^> ( 14 ) 
a 

where g a denotes the dimensionless conductance of the 
crth channel. 



For terms proportional to (N c \^) m g n with m 7^ n, the 
situation is quite different. Consider the N^g 2 term in 
Eq. (|l^). If we had suspended the sum over channels 
until the end of our calculation, we would have seen that 
these terms are proportional to 

[.9 2 ]tot = E-9*' ( 15 ) 

a 

Only when symmetry considerations constrain all the in- 
dividual channel conductances to be equal can we safely 
use [g 2 }tot = (gtot) 2 /N ch . 

Consequently, for the general situation in which the 
conductances in the separate channels are not necessar- 
ily equal, the appropriate equation for the fractional peak 
splitting is the following: 

/ a 0.1405( 5tot ) + 0.1491[.g 2 ] tot - 0.009798(. 9tot ) 2 + . . . 

(16) 

If we extended the expression to nth order in the di- 
mensionless conductances, it would contain factors such 
as 

[3 m ]tot = £.e, 

a 

where m < n and these factors appear both alone and 
in combination up to nth order in dimensionless conduc- 
tance. The details of the "fine structure" are fully char- 
acterized by the set of [(7™] tot for 1 < to < A c h, and 
the fractional peak splitting can be expressed in terms 
of these. Further modifications might be thought neces- 
sary to account for the "hyperfine structure" that results 
from allowing the tunneling amplitude t in Eq. (pj) to be 
a nontrivial function of ki and k2 . However, as long as 
the tunneling amplitude varies little over an energy range 
of order U2 , one would not expect Eq. (|l6|) to be changed 
substantially. 



III. THE STRONG-COUPLING LIMIT FOR 

iVcH =2 

The g 2 correction to the two-channel solution boosts 
confidence in the small-g end of our /-versus-g interpola- 
tion (see Fig. 3) but does little to improve the precision of 
theoretical predictions in the strong-coupling limit, a fact 
of particular concern for the experimentally relevant case 
of two interdot tunneling channels .ouliO The sections of 
the paper that follow improve the strong-coupling theory 
for N c b — 2 in two substantial ways. The first contribu- 
tion, presented here in Sec. Ill, is the calculation of the 
second term in the (1 — g) expansion about the g = 1 
ground state. This term, which is linear in (1 — g), is of 
interest both because it is significant in determining the 
shape of the /-versus-^ curve and because, in the cak 
culation that yields the primary (1 — g) ln(l — g) term,B 



G 



the (1 — g) term is naively ultraviolet-divergent. The sec- 
ond important contribution, which comes in Sec. IV, is 
the provision of powerful evidence that the coefficients 
of the leading terms in the (1—g) expansion are indeed 
independent of the high-energy structure of the theory. 

To calculate in the limit of g — > 1, we model the tun- 
neling link between the dots as a one-dimensional chan- 
nel with a delta function scattering potential at its cen- 
ter. This model was originally developed foij-lhe prob- 
lem of a single dot connected to a bulk leadBliil but was 
shown in Refs. 4 and 5 to be easily adaptable to that of 
a pair of coupled dots. Within this ansatz, the value of 
the double-dot charging energy is a simple reflection of 
the total number of electrons that have been transferred 
through this channel from one side of the barrier (dot 1) 
to the other (dot 2). In addition, as the system is effec- 
tively one-dimensional, the fermionic degrees of freedom 
can be bosonized, and the Euclideap_ action assumes a 
characteristic Luttinger-liquid formHOt 



S = 
Sq = 

Sint = 

5 b = 



Sq + Sint + Sb , 

^EEi w ™ii^)i 2 



ft 



dr 



VW 



rP 

V / dr cos[20F0 ct (t)]. 
_ Jo 



(17) 



In these formulas, 6 a (r) is a bosonic field that tracks the 
displacement of the one-dimensional electron gas at the 
barrier (x — 0), and CT (w m ) is its Fourier transform: 



r(r) 



ft 



r(w m ) , 



(18) 



where r is an imaginary time divided by H, ft is the inverse 
temperature (ft = 1/keT), and uj m is h times a bosonic 
Matsubara frequency (u> m = 2nm/j3). In addition, V is 
a measure of the barrier strength defined by V = for 
the delta- function potential VqS(x). As for the remaining 
parameters, vp is the Fermi velocity and, as in the weak- 
coupling theory, W is the bandwidth — the difference be- 
tween the lowest and highest single-particle energies in 
the channel. The inverse temperature ft will be taken to 
infinity in calculating the energy of the ground state. 

As in the weak-coupling theory, we ultimately want to 
parametrize the coupling between the dots by the dimen- 
sionless channel conductance g, rather than the barrier 
strength V. Accordingly, we need to find the relation 
between g and V. In our single-mode channel, g equals 
the single-particle transmission probability T, and (1—g) 
equals the reflection probability R. The leading depen- 
dence of the channel conductance on V equals what one 
would guess from the reflection probability of a single 
particle ttwiident upon a one-dimensional delta-function 
potentiaM 



(1 - g) = V 2 + 0(V 4 ) . 
Inverting this formula, we find that 

V 2 ^(l-g) + 0[(l-g) 2 } 



(19) 



(20) 



To lowest order, we have the approximation of Matveev,B 
V = a/1 — Sj which — it will be seen — is all that is re- 
quired for the calculations in this paper. 

Having prepared ourselves to switch from a solution 
in terms of V to one in terms of g, we proceed with 
the calculation of the ground-state energy. Our first 
move is to reorganize the action, expressing it in terms of 
bosonic fields that characterize the net charge and pseu- 
dospin degrees of freedom, where the pseudospin degrees 
of freedom correspond to "true spin" only if the two chan- 
nels correspond to spin-up and spin-down, respectively. 
Defining the charge field by 9 C = G\ + 62 + and 
the pseudospin field by 6$ = Q\ — 8 2 , we find 



S 



Sq ' + Sq ' + 5b • 



w m )| 2 , 



.--(c) 



^E(w + vO ^ c( " m)|2 ' 



vw 







dr cos 



V^0c(t) 
[V^0s(t) 



7T/9 



(21) 



The Euclidean action has now been written in terms 
of "high-energy" charge modes and "low-energy" pseu- 
dospin modes. We advance by integrating out the "high- 
energy" charge degrees of freedom. This integration is 
analogous to a renormalization in which one integrates 
out the highexjenergy degrees of freedom within a partic- 
ular channelO One begins with the generating functional 
for the Euclidean action of Eq. ( f2l| ) : 



Z = Jd[9 s (t)}Jd[6 c (i 



-S[e B (r),e c (r)] 



(22) 



where, as usual, time-ordering is implicit within the func- 
tional integral approach. One then performs the integra- 
tion over the fast modes to obtain the generating func- 
tional for an effective action depending only on the slow 
modes: 



D[9 s (t)} e 



-S e «[0 s (T)] 



s(T)] _ 



9 ° jD[e c ( T )]e- s o e- s * 
jD[0 c (T)]e-^ c) 



(23) 



Eq. (|2^) determines the effective action S e g. To solve 
for it, one Taylor-expands the exponential factor e~ Sb , 
performs the integral over charge degrees of freedom, and 
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re-exponentiates the result. Before doing any of this, it 
is useful to make the following definition: 



jD[e c ( T )]Ae' s o 



jD[6 c (r)]e-^ } 
One can then rewrite Eq. (^) as follows: 



(24) 



eft _ p~°0 



l-{S h ) c + -{S 2 ) c + 0(V 3 ) 



Upon re-exponentiation, one obtains 



(25) 



S cS = S ( S) + (S h ) c - \ ([5 b - (5 b )J 2 ) c + G(V 3 ) . (26) 



It is clear that to solve for the effective action to order 
V 2 = (1 — g), we must solve for both corrections to Sq 
on the right side of Eq. (p(J) . 

Details of the calculation of these terms are presented 
in Appendix B. The result is that 



where 



VW 



Sea — S, 



5, 



(i) 



S 



(2) 
b ' 



(|)| ? rfrco S [M(r)] , 



s^ = 



VW 



r*(°) / dn / dr 2 

Jo Jo 



x W(f ) 



I _ g— t#c(ti-t 2 ) 



- sinh[7rivr c (ri - r 2 )]} 
cos[Vtt6's(ti)] cos[V7r^ s (T 2 )] 



(27) 



The function JsTc(t) that appears in these formulas 
is the charge-channel correlation function, K c (t) — 
(0 c (t)9 c (0)) c - Its numerical value can be found from the 
formula 



K c {r) 



1. 



-Re 



duj ■ 



-(w+ lT ) w 



(28) 



To progress further, we define a new "unperturbed ac- 
tion" 5now = Sq +S^ . We then write down the Hamil- 
tonian that corresponds to this action: 



(i) 



I! TU b , 

' cos. 
7T V 2 



r SKo(°) cos (IP 



cos 



[^ B (0)] . (29) 



This is the Hamiltonian diagonalized by Matveev in 
Ref. 9 through a process of "debosonization" (see Ap- 
pendix B) in which the Hamiltonian is rewritten in terms 
of fermion operators dk and d: 



H - 



dk £, k d\d k , 



dk 



dl(d + dt) + (d + dt)d k 



(30) 



Here the single-particle energy the fermion interac- 
tion parameter A, and the wave- vector cut-off A have the 
formulas = hvpk, A = V cos(n p / 2) 2e~* hv p L^/tt 3 , 
and A = W/2hv F . 

Since the Hamiltonian is now quadratic in fermion op- 
erators, a Bogoliubov transformation brings it to the de- 
sired diagonal form: 



E 



(0) 
Now 



dk 



& (c[c k 



c\c k 



(31) 



where, if we write down only the terms of lowest order in 
V, replacing all others by an ellipsis, 



Cfe 
Cfc 



{d k + dl k ) 



1 

-^=(d k -dl k ) + 



(32) 



The correction to the V = ground-state energy is 
produced by the omitted terms in C k (for details, see 
Appendix B). In particular, using to represent the 
difference between the ground-state energies of H^ cw for 
V = and for arbitrary V, respectively, one finds that 



a£00 = 



4eyUo - 



2f>2 2( 1T P 
-V z COS ( -r- 



[ In 

-1 



V 2 



V 2 cos 2 



In 



(?) 

8e 7 ~ 



-In 



(33) 



As before, ip = W/U2, where the bandwidth W ^> U%. 

We see that the result for this first correction contains 
terms that are quadratic in V and logarithmically diver- 
gent in ip. This ultraviolet divergence is circumvented 
in Ref. 9 by the statement that one should replace W 
by U2 because keeping only the first term from charge- 
channel integration is only a good approximation for en- 
ergies less than the charging energy L^- The terms in 
Eq. ([33j) that are merely quadratic in V are thereby ren- 
dered finite and can be dropped in favor of the leading 
V 2 cos 2 (7rp/2) ln[V 2 cos 2 {-Kp/2)] dependence. 

To eliminate the logarithmic divergence more formally, 
one must calculate the shift in the ground-state energy 
that is induced by S<£' [recall Eq. As this term 

is itself quadratic in V and as we are only interested in 
knowing the ground-state energy to order V 2 , we can 

(2) 

drop all but the leading part of the -induced shift. 

(2) 

In expressing S^ in terms of the diagonalizing oper- 
ators of i?New) one may use the truncated formulas of 
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Eq. (32). The relevant shift in the ground-state energy is 

(2) 

then found by calculating the expectation value of 
in the ground state of Hj^ ew (see Appendix B): 



/ w 1 + 0.919(1 - g) ln(l - g) - 0.425(1 - g) 



(39) 



/o 



V 2 

l_ e -**.(#) 



(34) 



where units have been chosen such that h = 1 and terms 
independent of p have been dropped since they are not 
relevant to evaluation of the fractional peak splitting /. 
It is not too hard to see that the factor [1 - e -^K D {2x/W)^ 
in the integrand makes for an ultraviolet cut-off of order 
ip = W/U2 (see Appendix B). It is even easier to see that 
(1 — e~ x ) provides an infrared cut-off of order 1. Thus, 
one can surmise that the leading term from the integral 
is ln('i/)/2), which is precisely what is needed to cancel 

the ultraviolet divergence in A^ . 

What remains is to calculate the rest of the integral in 
Eq. EmI, which we call <E>: 



$ = lim 

Tp — >00 



dx 



In 



l_ e -^c(#) 



1 



(35) 



Numerical approximation of the integral in the limit 
V> -> 00 gives $ = 0.1703 ± 0.0002. 



One can now sum and A^ to get the strong- 

coupling energy shift through order V 2 . Having dropped 
terms that are independent of p, one has 



(2) 



A str (p) = ^t>W 



-1 + $ + ln 



(f ) 0° 



y 2 cos 



(36) 



We can now straightforwardly compute the fractional 
peak-splitting / in terms of the dimensionless conduc- 
tance g. As mentioned before, if we are only interested 
in obtaining the ground-state energy to order (1 — g), 
only the leading term of Eq. ( |20| ) is relevant in converting 
Eq. ( |36| ) to an expression in terms of (1 — g). The value 
of / = f P =i follows from the fact that, in the strong- 
coupling limit, 



f P = l 

In particular, Eqs. 

16e 7 



/ = 1 + 



7T 3 
7T ! 



Since $ » 0.1703, we have 



A str (p)- A str (0) 
C/ 2 p 2 /4 

and @ yield 



(1-9) + 



(37) 



(l- fl )ln(l-s) 







MS) 


- $ 







(38) 



Having determined the first corrections to the leading 
behaviors for both g — > and g — > 1, we now have a more 
plausible picture for the connection between the iV c h = 2 
weak- and strong-coupling limits (see Fig. 3). The fit 
to the data could be improved if the interdot ,-Cjapaci- 
tance were larger than experimentally estimated or if 
asymmetry between the dots were important.EfB In any 
case, whether or not such further emendations should be 
made, the theory is within the range of present experi- 
mental error. The corrections introduced in this paper 
have moved the weak- and strong-coupling predictions by 
reasonable amounts in the desired directions, increasing 
both the ease and the precision of interpolation between 
the weak- and strong-coupling limits. 
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Dimensionless interdot conductance g 

FIG. 3. Graph of the fractional Coulomb blockade con- 
ductance peak splitting / as a function of the dimen- 
sionless conductance per channel g in the weak- and 
strong-tunneling limits for N c h = 2. The new theoretical 
curves are depicted as solid lines. The old theoretical curves 
from Refs. 4 and 5 are dot-dashed lines. The dashed curve 
shows a possible interpolating function. Data points from 
Refs. 2 and 3 are given as triangles or stars; the two dif- 
ferent symbols correspond to different data sets. The value 
of / for the experimental data has been extracted from 
the measured splitting fraction /' by using the method dis- 
cussed in Ref. 4 with experimentally estimated values of 
20 aF for the constant interdot capacitance and 0.4 fF for 
the total single-dot capacitance. 2,3 



IV. INSENSITIVITY TO THE HIGH-ENERGY 
DENSITY OF STATES 

A. Insensitivity to Functional Form of Bosonic 
Cut-Off 

To have confidence that our coupled-dot calculations 
can be usefully compared to empirical data, we should 
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make sure that the result, / expressed as a function of 
g, is independent of the details of the band structure far 
from the Fermi surface, where the assumption of a con- 
stant density of states becomes invalid. We have done 
much to confirm such robustness in the regime of weak 
coupling, for we have shown there that, through second 
order in g, f(g) is independent of the bandwidth W and 
the filling fraction F as long as both FW and (1 — F)W 
are much larger than the charging energy 1/2- Such dual 
invariance indicates that we can simply shear off a non- 
trivial number of high-energy states without affecting 
the result. We would expect then that we could make 
less Draconian modifications of the high-energy density 
of states with similarly perfect impunity. 

With regard to the strong-coupling theory, matters 
have been left less assured. In Ref. 4, we introduced 
a factor of C2 multiplying the first term in Eq. (|3^) to 
guard against the possibility that the coefficient of the 
energy shift calculated in the Luttinger-liquid approach 
was partly a product of the approach itself and, in par- 
ticular, the manner in which the ultraviolet cut-off was 
imposed. Concern about such a possibility arises from 
the fact that the leading term in the (1 — g) expansion 
is proportional to the product of e~^ K "^ and A, where 
A is the generalization to no n-exp onential cut-offs of the 
normalization factor in Eq. ( B10| ) that gives the propor- 
tionality between the fermionic position operators and 
the exponentials of bosonic fields: 



^}(0,r) =Ae l ^f^ 



(40) 



Changing the nature of the bosonic cut-off [e.g., from the 
exponential e - "'"' to the Gaussian eT^I^ w ] causes 
the value of e~ 1 k K ^) to be multiplied by a constant fac- 
tor. Although one would hope that a similar shift in the 
value of A compensates for the change in e~^ Kc ^°\ to 
the authors knowledge, such a happy circumstance has 
not previously been checked to be true. 

Similar questions could be asked about the prefactor 
for the term linear in (1 — g), with which we associate 
a factor C3, where C3 = 1 for the Luttinger-liquid ap- 
proach with the standard exponential cut-off. This term 
is proportional both to |„4| 2 e~ 7r - ftrc (°) and to an integral 
that depe nds upon e~^ Kc ^ [see Eq. ( p4| ) in Sec. Ill and 
Eq. ( B26| ) in Appendix B]. Hence, in order to prove that 
the two leading strong-coupling terms do not vary with 
the choice of cut-off function, one must show that neither 



the product V x = \Afe 



-7rK c (0) 



nor the integral 



V, 



1 



j_ r pw ' 2 (f3W 



dx 



\ 2 



1 



-wK c (2x/W) 



(41) 



assumes different values when the shape of the cut-off is 
changed. Though we do not have a general proof that T>\ 
and T>2 are independent of the cut-off function, we can 



show that they remain the same for a whole class of func- 
tions that includes the exponential cut-off and that they 
are similarly unchanged when one replaces the exponen- 
tial cut-off by a Gaussian. We believe that these facts 
are convincing evidence that the prefactors in Eq. ( |38| ) 
are insensitive to the nature of the high-energy cut-off. 

First, we prove that T>\ and T>2 are the same for all 
cut-offs of the form 



M 



(42) 



where either M is finite or, for large m, b m falls to zero 
faster than m _l »/m! for some real C > 0. As usual, it is 
assumed that aU2 <C 1, where a — 2/W . We add the 
further assumption that [(to — 1)! b rn aU2] <C 1 for all to. 

The first step in our proof is to solve for the change 
in g-^^^o) wneri one goes from the standard expo- 
nential cut-off v{oj, a, {0}) to the more general form 
v(ui,a,{b m }). The formula for K c (t) [recall Eq. (J2S* 
becomes 



1 f°° 

K c (t) = -Re / dio v(u),a, {b m }) 

T JO 

We can write the change in K c (0) as 

M 



2U 2 



sk c (o) = i y> 

TV * 



duj- 
w 



2U 2 



(43) 



(44) 



Using to = (uj + 21/2/^ — 2U~2/tt) and the binomial theo- 
rem, we can expand uj m in powers of (to + 21/2/^)- The 
integration is then straightforward and yields 



M 

E 

m— 1 



b m (m-l)l[l + 0(aU 2 )]. 



Dropping the correction, we have the result 



-J2b m (m-iy. 



-7i-if c ,o(0) 



(45) 



(46) 



where K c q (t) is the correlation function for the standard 
exponential cut-off. 

Calculation of the change in the normalization cop-, 
stant A is more complicated. Following V. J. EmeryEHl 
we find 



dx ^^{b m }) ^' h J F 1 + c c ^ 



A bit of calculation reveals that 



\A\- 



dx 



a — ix 

£ b m (m-iy.(^) 



It is apparent that the bracketed factor exa ctly cancels 
the factor that multiplies e ~ 7rK "- 0< - ' ) in Eq. (||). Thus, 
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in order for T>\ = \A\ 2 e~ 7 ' Kc ^ to be unaltered, the value 
of the integral in Eq. (fi§| ) cannot change as the b m are 
varied. In short, the partial derivative of the integral 
with respect to each of these coefficients must be zero. 
The partial derivative with respect to b m is given by the 
following formula: 



P m =(m-1)\ [°° dx(-^-r-\ 



x e 



a — ix / 
EM"»-i)«(3^) T 



m+l 



+ C.C. 



(49) 



Let z — a I [a — ix). The resulting integral in the com- 
plex z plane follows a closed path, beginning and ending 
at z = 0: 



ia(m-l)! <bdzz m - 1 e^ bm{m - 1) - zm + c.c. (50) 



P„ 



For M finite or b m falling off faster than m~^/m!, the 
integrand is analytic throughout the region enclosed by 
the contour. By Cauchy's Theorem, P m = 0. 

Since the integral of Eq. (|8|) does not vary with b m , 
we can make the statement 



\Af 



-1)! 



\A 



(51) 



where Ao is the normalization factor for the standard 
exponential cut-off. Eqs. (46) and ( pl| ) yield 



V 1 = \A\ 2 e-* K ^ = | A |2 e 



-77if c , (0) 



(52) 



We have now shown that, for the class of cut-offs 
v(uj,a, {& m }), C2 is constant. 

What about C3? To determine its fate, we must find 
the change in the quantity T>2 [recall Eq. (f4l|)]. After 
substituting a for (2/W), we follow essentially the same 
path that we blazed in determining the change in K c (0) 
and find that the change in K c (ax) is given by the for- 
mula 



1 



M 



5K c (ax) = -j- ^ b m {m - 1)! 

rn — 1 



1 



(i 



+ c.c. 



■ IX) 



(53) 



Employing this, we can break the integral on the right 
side of Eq. ( f4"l| ) into two parts, the first of which is from 
to V 1_e , where < e < 1 and V = (W/U 2 ) > 1. In 
this interval, the contribution from the entire term pro- 
portional to e -^ K ^( ax ) can be shown to be zero in the 
limit ip — > 00. In the remainder, 5K c (ax) is on the or- 
der of 1/x 2 , which implies that the correction due to the 
generalization of v(uj, a, {b m }) is proportional to 



/3W/2 
j/i 1 " 



dx < 



pw ' 2 dx 



l/l 1 



which also equals zero in the limit tj) — > 00. Therefore, T>2 
is constant, and we have proven that our strong-coupling 



results are insensitive to varying the cut-offs within the 
class a, {b m }). 

The values of D\ and T>2 can be shown to be similarly 
unaltered when we switch from the exponential cut-off to 
a Gaussian: 



v g (uj, a) 



e 4' 



(54) 



Solving for e 7rKc (°) with this weight function, one dis- 
covers that 



(55) 



where 7 is once again the Euler-Mascheroni constant. 
The normalization coefficient Ag has not been solved for 



analytically. However, starting from Eq. (47), one finds 
that 



\A C 



tux I dxcos 
'0 

x e 



-Erf(x/2) 



f*dye~y 2 ^Er{(iy/2) 



(56) 



where Erf(x) = (2/^/tt) fydte~ t is the error function. It 
has been confirmed numerically that through at least 12 
digits the product T>± t a = \AG\ 2 e~* Kc ' G ^ agrees with 
the exponential cut-off. By arguments similar to those 
used for the class of cut-off functions studied above, it 
has also been shown that in the limit tp — » 00, the integral 
T>2.g is the same as for the exponential. The coefficients 
in Eq. ( p8| ) are again unaltered, and it seems reasonable 
to suppose that the invariance is general. 



B. Insensitivity to Fermionic Filling Fraction 

Thus, it appears fairly certain that modifying the high- 
energy density of states in the bosonized theory does 
not affect the results of Sec. III. Nevertheless, having 
solved the weak-coupling model for the general case of a 
fermionic system not necessarily at half-filling and having 
seen that, when expressed as functions of the tunneling 
amplitude, both the conductance and the fractional peak 
splitting depend upon the filling fraction, one might won- 
der what happens to the strong-coupling results when 
one begins with a fermionic system that is not neces- 
sarily at half-filling. Since Luttinger-style bosonization 
assumes symmetry between occupied and empty states, 
such a system can only be properly bosonized after the 
asymmetric fermion states have been integrated out. For 
example, if the system is below half-filling [F < (1 — F)] 
and the zero of energy is at the Fermi surface (ep = 0), 
the fermionic single-particle states with energies between 
FW and (1 — F)W must be integrated out, leaving a sym- 
metric effective theory with single-particle energies rang- 
ing from —FW to FW. Only after this symmetrization 
can the theory be bosonized without losing knowledge of 
the fermionic filling fraction F. 
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The task before us, therefore, is to "symmetrize" the 
fermionic theory that lies behind the bosonized action of 
Eq. (0) . The archetypal fermionic Hamiltonian consists 
of the usual three parts: the single-particle kinetic ener- 
gies, the multiparticle potential energy, and the backscat- 
tcring barrier. The Hamiltonian therefore takes the fol- 
lowing form: 



H 



Hr + He + Hb , 



3=1 a k 

H c = U 2 (n - p/2) 2 , 

#b = £ £ w (4k 2 <T C iki<r + H.c.) , 



(57) 



where = twpk and j is the index that distinguishes 
between right-movers (j = 1) and left-movers (J = 2). 

The operator h is now somewhat more complicated 
than in the weak-coupling theory. In its simplest form, 
it can be written as 



2 ^ 

j=i 



dx [8 (at) - Q(-x)] ipUx^jix) 



(58) 



where Q(x) is the Heaviside step function and t/jj is the 
annihilation operator in position space for a right-moving 
(J = 1) or a left-moving (j = 2) fermion. After writing 
the components of the integrand in the momentum rep- 
resentation and integrating over x, one hnds that, for a 
one-dimensional system of length L, 

"EEE^d-w, (59) 



L ^ ^ ^ k 2 

j <r k\ki 



ki 



which is_aquivalent to the integral version obtained by 
MatveevEJ from the observation that dh/dt equals the 
current operator at x = 0, the point of "division" be- 
tween the two dots. (This point is, of course, not entirely 
well-defined in the limit g — > 1.) 

The above equations for the Hamiltonian and number 
operator are presented as discrete sums. For future refer- 
ence in implementing the symmetrization of the theory, 
we write the components of our Hamiltonian in integral 
form: 



H 



K 



2 

£ 

H c = U 2 (n-p/2) 2 , 



( Hvf 



dk^k Cjka C jk<J ) 



H B =v (^^j J dk i Jdk 2 (4 k2a c lkia + H.c.) (60) 

where S is the level-spacing for the one-dimensional sys- 
tem (6 = 2nhvF j ' L) and 



c jkic 



k 2 - h 



(61) 



For the fermionic strong-coupling model of Eq. Q57^ 
calculation of the channel conductance between the dots 
proceeds along the same lines as for weak coupling (see 
Ref. 4). In fact, since the density of states is con- 
stant in both theories, setting U 2 — — the first step 
in the conductance calculation — renders them essentially 
identical, the only differences being in the last term, 
where t has been replaced by v and the index i for dot- 
1 or dot-2 fermions has been replaced by the index j 
for right-movers or left-movers. Accordingly, unlike the 
weak-tunneling term H?, the perturbation Hg scatters 
fermions backward instead of transporting them forward 
and therefore causes a reduction in the conductance of 
the unperturbed system. Recalling the size of the con- 
ductance induced by Ht in the weak-coupling model, it 
is not hard to see that the channel conductance in the 
strong-coupling model is given by 



4x 



\l + (l + i7 1 fxV 



(62) 



where x — (irv/5) 2 and 77 = (1/tt) \n[F/ (1—F)]. As in the 
weak-coupling theory, the result becomes troublesome as 
X becomes large. However, we should be able to trust 
its testimony that the filling fraction does not affect the 
interdot conductance through second-order in (ttv/S). 

This is all we need to know, for (ttv/S) can be 
straightforwardly written in terms of our previous strong- 
coupling parameter V. The relation is V = 2(irv/8), and 
it follows that x = (V /2) 2 . We recover the leading-order 
result of Eq. (|20| ) and see that, to order V 2 , the channel 
conductance g is independent of the fermionic filling frac- 
tion F. If we can likewise show that the relation between 
V and the differential energy shift 



5A(p) = [A stI (p) - A str (0)] 



(63) 



does not depend on the filling fraction F, we will know 
that the same is true for our final strong-coupling result, 
the expression for f(g) in Eq. d38|). 

To prove SA'a invariance with respect to F, we sym- 
metrize the fermionic theory through a renormalization 
in which we integrate out all single-particle states at an 
energy distance of W/2 or more from the Fermi sur- 
face, where U 2 -C W <C W . The resulting symmet- 
ric theory with bandwidth W' can be bosonized without 
further qualm. However, as renormalization generates 
terms that are not present in the original Hamiltonian, 
we must check to see what relevant effects these have 
upon the low-energy theory. We must also keep track of 
any contributions to i5A that arise from the high-energy 
degrees of freedom alone. 

Before we go about doing this, a comment on our ap- 
proach is in order. One might view the proposed renor- 
malization as occurring in two distinct stages: first, we 
integrate out the asymmetric particle-hole states; then, 
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we integrate both particle and hole states down to en- 
ergy W' . Since all that we will need to consider are the 
general scaling properties of the terms generated during 
the renormalization process, the distinction between the 
stages is of no importance and is henceforth ignored. 

The argument resumes. Since our interest is in the 
Coulomb blockade, the renormalization scheme we use is 
designed to leave the Coulombic interaction term Hq un- 
changed. After wave-vectors between the original wave- 
vector cut-off A and the new wave- vector cut-off A/6 
(where b > 1) have been integrated out, the theory is 
re-scaled by writing it in terms of a new set of wave- 
vectors ki, — bk. fnvariance of He is achieved by re- 
scaling the fermion creation and annihilation operators 



as well: c 



t 

jk h a 



-VM. 



-jka ■ 



[One might prefer to say that 



the coherent-state Grassman variables that correspond 
to the operators are re-scaled (see Ref. 38).] The effect 
of renormalization upon the parameters {Kvf 1 5), U2, and 
v(Hvf/5) of Eq. (|5^) is as follows: 





1 

= b- 1 


hvp 


s 




s 



[u 2 ]' = u 2 , 

[v(hv F /5)}' = b- 1 [v{hv F /8)\ 



(64) 



The backscattering Hb is revealed to be dangerously ir- 
relevant. Though it scales like an irrelevant term, we 
cannot safely set it to zero as we know from Eq. (|3^ ) 
that the energy shift is singular as v — > 0. 

In addition to rescaling the terms in the original Hamil- 
tonian, renormalization generates terms of its own. It is 
not hard to see that all but the new backscattering terms 
are irrelevant. The original Hamiltonian consists of the 
kinetic energy Hk , a two-body interaction Hc%, a one- 
body interaction Hqi, and a backscattering term Hb- 
Hc2 and Hex are normal-ordered operators given by the 
following formulas: 



Hc2 — 



31 ,J2 o\,a% 



Hex — 



c jk 2 <J C jki a ' 



(foj - h)(k2 - kx) 

A-fc 2 



x In 



3 o 

A-fci 



A + fci 



- In 



A + fc 2 



, (65) 



where A = W/Hvf- Hc2, Hex, and Hb can be repre- 
sented by Feynman graphs [see Fig. 4(a)], which can then 
be connected to construct the terms that renormalization 
adds to the Hamiltonian. As usual, the internal lines of 
the second-generation graphs carry only high-energy mo- 
menta which lie within the shell of wave- vectors that are 
integrated out. 





(3) 



(3) 



(1) 





(C) 



(2) 



(d) 



FIG. 4. Feynman diagrams for integrating out sin- 
gle-particle energies a distance greater than W' from 
the Fermi surface in the fermionic version of the 
strong-coupling theory, (a) The three building-block Feyn- 
man graphs. Diagram 1 corresponds to the two-body 
Coulomb interaction Hci- Diagram 2 corresponds to 
the one-body Coulomb term Hoi- Diagram 3 represents 
the backscattering Hb- (b) Second-generation m-body 
graphs constructed by contracting Hci 's and Hci 's. These 
terms are all irrelevant to the low-energy theory, scaling to 
zero under renormalization. (c) Second-generation graphs 
formed from combinations of Hb, Hc2, and Hci- Terms 
such as Diagram 1 that contain an even number of Hb's are 
irrelevant under renormalization. Diagrams 2 and 3 involve 
odd numbers of Hb's and are therefore dangerously irrel- 
evant. Nevertheless, they are negligible in size compared 
to corresponding low-energy graphs and therefore can be 
safely discarded, (d) An example of a closed diagram used 
to calculate the contribution to the energy shift from the 
degrees of freedom that correspond to single-particle ener- 
gies more than W' from the Fermi surface. 



Given such rules for making second-generation terms, 
one can deduce that, whenever one creates a new term 
by connecting lines emanating from the He 2 and Hex 
graphs [see Fig. 4(b) for examples], one picks up a scaling 
factor of For example, Diagram 2 of Fig. 4(b) repre- 
sents a two-body interaction produced by contracting one 
Hci with one Hex- This new interaction term is similar 
to Hci except that the denominator contains only one 
power of (fej — £3) or (& 2 — &i) and, consequently, is less 
singular than Hc 2l which is fixed under rescaling. Thus, 
the second-generation term must shrink under renormal- 
ization. Indeed, all such graphs formed from contracting 
the Coulombic interaction terms are similarly irrelevant 
and scale to zero under renormalization. They can be 
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ignored in the effective theory. We should expect this 
result. Otherwise, our Coulomb blockade model would 
probably never have been useful at all. 

As for graphs that involve the backscattering term Hb 
[see Fig. 4(c)], we need only consider these to order 
v 2 , for we go no further in calculating /(<?). Depending 
upon how many Coulombic interaction terms are intro- 
duced, the second-generation graphs that contain Hb all 
scale down by at least a factor of b^ 1 . Consequently, 
all but those which contribute to low-energy backscat- 
tering are irrelevant. Thus, we can drop graphs such as 
Diagram 1 of Fig. 4(c) that contain an even number of 
Hb's. Graphs containing an odd number of Hb's are 
dangerously irrelevant but can ultimately be ignored be- 
cause they are negligible compared to the corresponding 
graphs that can be constructed from the low-energy por- 
tions of the original Hc2, Hex, and Hb- Diagrams 2 and 
3 of Fig. 4(c), for example, are of order v(U2/W). If we 
had renormalized down to W" , where Ui <§C W" <C W, 
we would have found the corresponding graphs to be of 
order v(U2/W"). The contribution from energies above 
W' is therefore seen to be merely perturbative in relation 
to the contribution from energies between W" and W' . 
The conclusion is that we can drop the parts of the graphs 
produced by integrating over energies greater than W. 
Returning to our original renormalization down to W', 
we see that the graphs produced here have been shown 
to be negligible. The argument that the symmetrizing 
renormalization does not cause any significant changes 
in the low-energy Hamiltonian is complete. 

Having disposed of the concern that the process of 
symmetrization might leave us with important new low- 
energy terms, we now show that any constant terms pro- 
duced are similarly insignificant. Such constant terms 
correspond to closed diagrams constructed from the orig- 
inal Feynman graphs. Since all lines are internal, they 
all represent the propagation of high-energy excitations. 
There are obviously an infinite number of closed dia- 
grams. Fortunately, we can limit our attention to a 
certain subset. We need not concern ourselves with di- 
agrams involving less than two Hex graphs: diagrams 
with only one Hex graph must sum to zero as A(p) is 
even in p; diagrams with zero Hci graphs cannot con- 
tribute to the differential energy shift 5A. Similarly, in 
any pertinent closed graph, Hb must appear a nonzero 
and even number of times. It cannot be absent as terms 
that do not include it shift all relevant ground-state ener- 
gies equally and are therefore unimportant. Furthermore, 
in any closed graph, it must appear an even number of 
times because Hb is the only term that exchanges right - 
and left-movers. Thus, all the diagrams we need con- 
sider consist of a nonzero and even number of Hb's, at 
least two Hqi's, and an arbitrary number of flea's [see 
Fig. 4(d) for a canonical example]. 

Each such diagram corresponds to a number of time- 
ordered terms in Rayleigh-Schrodinger perturbation the- 
ory. For a Feynman diagram with r internal lines, the as- 
sociated Rayleigh-Schrodinger terms have r integrations 



over momenta and (r — 1) propagators with denomina- 
tors linear in the momenta. If the Hq2 graph appears 
m,2 times in the Feynman diagram and the Hci graph 
appears mi times, there are m = (2m,2 + mi) > 2 addi- 
tional denominators linear in the momenta, which have 
their origin in the wave-vector denominator of h [recall 
Eqs. © and ©]. The propagator denominators are al- 
ways on the order of W or greater. The fi denominators 
are of the form (k — k'), where k and k' are both in the 
high-energy wave-vector shell. Thus, these denominators 
can go to zero. However, the contribution from the re- 
gions where they become zero is negligible, the somewhat 
simplified explanation being that, when one of them goes 
to zero, the rest of the integrand can be treated as essen- 
tially constant, and we have 

where A' = W'/hvp and the constant e -C 1. It follows 
that contributions to the overall result only come when 
the h denominators are themselves of order W /Hvp. 

As a result, what remains is a nonsingular integration 
over r momenta of an integrand that is proportional to 
[hi . . . kr r -i+ m )]~ , where the ki are linear in the mo- 
menta over which we integrate. Noting that the only 
other momentum dependence comes from the logarith- 
mic term of Hci, we see that, in energy units, the result 
of the integration is of order (l/W 1 )" 1-1 . We now multi- 
ply the result of our integration by the various factors of 
U2, v, and 5 that stand aside the integral. For a closed di- 
agram in which Hb, Hc2, and Hex appear j, m%, and mi 
times, respectively, the contribution to the energy shift is 
readily seen to be of order Uziv j '<5) 2j {U2 /W') m ~ 1 , where 
m = (2m 2 + mi) > 2. As (U 2 /W) < 1 and the overall 
energy shift is of order U2, these terms are negligible. 

Thus, at least to order v 2 , integrating out all particle 
and hole excitations at distances greater than W' /2 from 
the Fermi surface produces neither relevant new terms 
in the low-energy Hamiltonian nor significant constant 
contributions to the differential energy shift. As what 
remains is a fermionic theory at half-filling, the result for 
f(g) in Sec. Ill is unaffected by possible "high-energy" 
deviations from this condition, an important property if 
we wish to compare our predictions with empirical data. 
We would hope that a similarly universal solution for 
f(g) could be found to higher orders in (1 — g). How- 
ever, if the formula for the interdot conductance [recall 
Eq. (|S2|) ] is correct to some non-leading order, such over- 
all independence of the filling fraction must — as in the 
weak-coupling limit — come through cancellation of the 
separate filling-fraction dependences of the conductance 
and the energy shift when one is expressed in terms of 
the other. If this were shown to be true, we would see 
once again that the interdot conductance g and not the 
bare matrix element for tunneling or reflection is the cor- 
rect parameter to achieve a universal description of the 
coupling dependence of a double-dot Coulomb blockade. 
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V. CONCLUSION 

The present paper substantially improves the results of 
earlier theoretical work on the CoulombJpleckade peak- 
splitting for two coupled quantum dotsoO making an 
important contribution to the growing body of theo- 
retical .and ^experimental work on such coupled-dot sys- 
temslflii3cj~E3. By extending the weak-coupling theory 
to second order in g for arbitrary iV c h, it has shown how 
the small- iV c h theory crosses over to the large- theory 
in the neighborhood of w 10. Furthermore, it has 
demonstrated that, at least for the leading two terms in 
the weak-coupling theory, the channel conductance g is 
the "correct" parameter to use in constructing a theory 
for the peak splitting that is universal in the sense that 
it is does not depend on the high-energy band structure. 
Finally, this paper has made the N c ^ = 2 theory both 
stronger and broader — broader in that the sub-leading 



term is calculated; stronger in that the leading and sub- 
leading terms for strong-coupling are confirmed to be in- 
sensitive to the manner in which the high-energy cut-off 
is taken. Thus, the paper has made more plausible ef- 
forts to connect weak- and strong-coupling behaviors and 
to compare theoretical resjilts,with the data from recent 
two-channel experimentsBfali2l 
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APPENDIX A: DETAILS OF THE WEAK-COUPLING CALCULATION 



As described in Sec. II, the procedure in evaluating the fourth-order energy shift is to calculate the (7V c h) 2 and N c h 
terms separately. Calculation of the (iV c h) 2 terms is facilitated by rewriting them in terms of two energy variables 
instead of four. Calculation of the iV c h terms is made easier by differentiating twice with respect to p while performing 
the integrations over energy and then integrating twice with respect to p at the end. Terms that are constant or linear 
with respect to p cancel in the final result, the relative energy shift (Ao — A p ), so we have not lost useful information 
as a result of the double differentiation. 

As mentioned in Sec. II, the "wrinkle" in these computations, the appearance of integrals of the form 



ax 



B) 



A 



is resolved by Taylor-expanding the logarithm about (x + A) for (B — A) < \x + A\ and about (B — A) for (B — A) > 
\x + A\. For A < 0, one first breaks the integral into the intervals (0, \A\ — e) and (\A\ + e,Rip). After this, one 
proceeds as usual. The results are 



V 



m , ln(:r 
ax 



.4 



L) =iln 2 (i# 



A) In 2 A 

' 2 



E 



(-i) 



n+l 



B-A 
A 

for < (B 



1 _ f B-A 
\ Rip + A / 

A) < A < (Rip + A), 



i ln 2 (i?V + A) - \n(B - A)lnA + ^ ln(B - A) 



E 



(-i) 



n+l 



A 



B-A 



/ B-A 
\Rf+A 



71" 



for < A < {B - A) < (Rip + A), 



= \n(B - A)[ln(Rip + A) - In A] 



E 



(-i) 



n+l 



Rlp + A 



A 



B-A J \B-A / 
for < A < (Rip + A) < (B - A), 
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i ln 2 (i# - \A\) - ln(B + |A|) In \A\ + ± ln 2 (B + + y 

^ 1 / H \ n ~ (-l)n+i / ^ + |A1 
^n 2 ls + |A|/ ?i 2 Ie^-UI 

n— 1 x ii/ n— 1 v 11 

for A < < (B + \A\) < (Rip - \A\), 



= \n(B + \A\)[ln(Rip - \A\) - In \A\] 



^n 2 VS+U|/ ^ n 2 \B+\A\ 

for A < < {Rip - \A\) <(B+ \A\). (Al) 

These five integrals are all we need. In confirming that-the solution for (B — A) < (Rip + A) evolves continuously into 
that for (Rip + A) < (B — A), it is useful to rccognizdI3 that 

n=l 

Having equipped ourselves to smooth the "wrinkles," we can proceed with a fuller description of calculation of the 
fourth-order terms. The (A c h) 2 calculation is reviewed first. An illustrative segment of the calculation follows. 

In Sec. II, it was remarked that each of the (A c h) 2 terms could be written in terms of two energy variables 
(ej = ea — ei, 64 — e 2 ) instead of the original four. The "cost" of this conversion is the appearance of a nontrivial 
density of states: 

r&f r^F pW^ p\fy r^-F r^F 

I dei / de 2 de 3 de± h(e 3 - ei, e 4 - e 2 ) = / dt\v(ei) I de n v(eu) h(e h e n ), (A3) 

JO JO J( F JtF Jo Jo 

where v(e) is the density of states. For a system at or below half-filling, 

{e for < e < e F , 

e F for e F < e < (W - e F ), (A4) 

(W-e) for (W- e F ) < e < W. 

(We need not worry about a system above half-filling as such as system can be mapped to one below half-filling 
through an exchange of particles and holes.) 

Using the new variables t\ and en, we can sum the integrands for the (iV c h) 2 terms shown in Fig. 2 (the others are 
obtained by taking p — > —p). If drop the common factor —N c \ l (t/5) 4 U2, the result is the following: 



[ei + U 2 (l - p)} 2 hi + U 2 (l - p)][ea + d + U 2 (A - 2p)} 
2 

+ fa + U 2 (l - p)f[e n + U 2 (l + p)][e u + ej] ' (A5) 



It is not hard to show that 



j dei / den [htotfa, e n, p) ~ /Hot fa) en, 0)] = 

Jep Jew 



1 e F Je F 

in the limit ip = W/U 2 — ► 00. Accordingly, we need only calculate 

/ dei deuh tot (ei,eu,p). 
Jo Jo 

The process of evaluating t his double integral is lengthy but straightforward. The only "wrinkles" that appear — 
integrals of the form of Eq. Al — are no longer problematic. The end result is Eq. (||) of Sec. II. 

We now move to consideration of the fourth-order terms linear in the number of conducting channels. Recall that 
the (A^h) 2 terms were added before the integrations over energy were performed. This order of tasks is reversed for 
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the iVch terms, the computation of which revolves primarily about finding a favorable permutation of the operations 
of differentiating and integrating with respect to p, integrating with respect to the ith energy variable, and integrating 
by parts. Consequently, perhaps the best way to describe the derivation of the N c h contribution is to walk through 
the computation of a single illustrative term. After seeing the methodology employed in calculating this term, the 
tireless reader should have little difficulty in computing the rest. 

The representative term we choose is that corresponding to Diagram 2 of Fig. 2(b). This term involves an exchange 
of a pair of electrons and, consequently, picks up an "exchange" minus sign. The diagram depicts the following 
sequence of events: 

I. Electron 1 tunnels from dot 1 to dot 2, going from a single-particle state with kinetic energy e\ to one 
with kinetic energy £3. The energy of the resulting double-dot state relative to that of the unperturbed 
ground state is [€3 — £1 + U 2 (l — p)]- 

II. Electron 2 tunnels from dot 1 to dot 2, going from a single-particle state with kinetic energy £2 to 
one with kinetic energy £4. The system's energy relative to the unperturbed ground state is now 
[e 4 + e3-e 2 -ei + 2[/2(2-p)]. 

III. Electron 2 tunnels back to dot 1, settling into the initial single-particle state of Electron 1. The ensuing 
relative system energy is [£3 — e 2 + Uz(l — p)]. 

IV. Electron 1 tunnels back to dot 1, settling into Electron 2's initial single-particle state. The unperturbed 
ground state has been recovered. 

With all the intermediate-state energies known, it is easy to write down the contribution to the fourth-order energy 
shift: 

4S» = * 4 £££ 



£3 - £2 + U 2 (l - p)\ 
a £1,62 £3,64 

x \ x (A6) 

[e4 + e 3 -e2-ei+2C/ 2 (2-/9)] [e 3 - ei + U 2 (l - p)} ' 

The sums over t\ and £2 extend from to the Fermi energy 6f ■ Those for £3 and £4 go from to the bandwidth 
W. The sum over the channel index a results from the fact that Electrons 1 and 2 can share any one of the N c b 
tunneling channels. Though the formula contains such unphysical terms as that for which t\ = e 2 , such terms are 
down by factors of the level spacing S divided by S/FW or 6/(1 — F)W, and their inclusion has no effect in the limit 
W/5 -> 00. 

Accordingly, we can cease worrying about these terms, for we a ssume that 8 <C U 2 <C W, a postulate that permits 
us to work in the continuum limit, replacing the sums in Eq. (A6) by integrals: 

/i\ 4 f eF f eF f w f w 1 

aZm=»*(- s ) 1 *.y„ *>i **i *■ [«-,, +„,(!-,)] 

x - x - (A7) 

[e* + ea - e 2 - ei + 2U 2 (2 - p)] [e 3 - ei + U 2 (l - p)] ' V ; 

These integrals can be rewritten in terms of dimensionless variables x$ : 

for i = 1 or 2 



u 3 



u 2 



^ for i = 3 or 4 . 



(A8) 



With this choice of integration variables, it becomes clear that A]y 2 (p) is linear in U 2 . Specifically, we find that 



^,2(p)=^ch(0 U 2 xI(p,F,i>) 



I(p,F,ip) — I dx\ I dx 2 / dxz / dx4 



[x 3 + x 2 + 1 - p] 

x -r x . (A9) 

[X4 + X3 + x 2 + xi + 2(2 — p)\ [X3 + xi + 1 — p\ 

All the shuffling of notatation still leaves us confronted with a quadruple integral. Opting to postpone a frontal 
assault, we try a sidestepping movement, computing the partial derivative with respect to p: 
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Fp 



dxi 



Fp 



dx- 



(l-F)V 



(l-F)i/> 



dxd 



][ 



][ 



][ ][ 



(AlO) 



where the subscript p signifies that I p is the partial derivative of / with respec t to p and the brackets on the right-hand 
side of the equation have the same contents in the same order as those in Eq. (A9). As the third term in the integrand 
of Eq. (AlO) differs from the first only by an exchange of the indices 1 and 2, we can drop the third term and double 
the first. When the enhanced first term is integrated by parts with respect to #2, the products are two triple-integral 
terms and a quadruple-integral term that exactly cancels the second term of Eq. (AlO). With the definitions Aq = 
and Ai — Ftp, we have 



p=0 



dx\ 



(l-F)i, 



dx- 



(l-F)i/> 



dxA 



[x 3 + A p + 1 - p] 
1 



[xa + x 3 + xi + A p + 2(2 - p)\ [x 3 + x x + 1 - p] 



(All) 



Having benefited once from differentiation with respect to p, we try it again. The second derivative of / with respect 
to p has the following form: 



J — 2 
1 pp ~ z 



Fiji 



dx\ 



(I-F)tP 



dx 3 



(1-F)V> 



dxA 



1 



1 



][ 



][ 



][ ][ 



(A12) 



where the bracket contents correspond — in order of appearance — to those of Eq. (All). I pp lacks the convenient 
symmetry between first and third terms that was so handy before. Nevertheless, integration of the first term by parts 
with respect to X3 still helps. The triple integrals that result cancel the third term and half the middle term, leaving 



2 EE 



(-1) 



p+q 



p=0 q=0 P 



A p + B q + 1 - p 



Fiji 



dx\ 



(\-F)i> 



dxA 



[xa, + xx + A p + B q + 2(2 - p)] [xi+B q + l- p] 



2V(-1) P / dxi / dx 3 / dx 4 

p =o J° J° -'o 



1 



[xa + A p + 1 - p] 



[xi + x 3 + xi + A p + 2(2 - p)] 2 [x 3 +xi + 1 - p] 



where Bq = and B± = (1 — F)ip. 

We now straightforwardly integrate over X4, using the relation 

1 1/1 



The result is that 



(x + a) (x + b) b — a \x + a x + b 



^ ^ ^ Ap + B„ + 1 - p Jn 

p=0q=0r=0 P 1 1 ' 



F ^ ln[x + A p + B q + B r + 2(2 - p)] 



x + B q + 1 - p 



p^UAp + ^ + Z-pJo 



Fij] 



dx\ 



(l-F)V 



dx 3 



[x 3 + A p + 1 - p] 



[xa+xi + l- p] [x 3 + xi + A p + B q + 2(2 - p)\ 



(A13) 



(A14) 



Recalling Eq. (Al), we see that, as if> — > 00, the leading part of the first term in Eq. ( |A14 ) goes as [\n 2 ^/(A p + B q + 
1 — p)] and therefore goes to zero unless p — q = 0. The same is true for the second term — which upon integration over 
x\ will have a form like that of the first term. Hence, we can eliminate the sums over p and q and, after integrating 
the second term over x±, have 
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_ 2 * p ln[s + g r + 2(2-p)] 



x+A r +l-p 



x+A,,+2(2-p) 



X + 1 - p 



(A15) 



We recognize that the second term is nontrivial only for r — and apply Eq. (Al) to do the remaining integrations 
over x. After dropping terms that go to zero as ip — > oo, we arrive at the "final" formula for I pp : 



I - + J (2) 
1 pp — 1 pp ' 2 pp 



2 I pp 



1 OO 

ln([l - F]ip) In(Ftp) - - ln 2 (F^) + ^ 



(-1) 



n+l 



F 



1 - F 



TV 



ln([l - F]1>) ln(l - p) - - ln 2 (3 - p) + ln(3 - p) ln(l - p) 



(-1) 



n+l 



n=l 



\_-p_ 

3-p, 
for F < (1 - F), 



^ OO 

-ln 2 ([l-F]^)-^ 



(_l)n+l f\ — F 



F 



for F > (1 - F); 



ln([l-F]V0 In(l-p) 



iln 2 (3-p) + ln(3-p)ln(l-p) + £ 



(-1) 



n+l 



1-P 

3-p 



3-p 



r(2) 
pp 



^ M 1 1 2 



1-p 

3-p 



E 



(-i) r 



3-p 



(A16) 



Before undoing the differentiations with respect to p, we pause to remark on the meaning that can be attached to the 
derivatives I pp and I p . The second derivative I pp can be interpreted physically [after multiplication by N c \JJ2(t/5) A ] 

as reflecting a change in the effective differential charging energy J7 c ff = 2[d 2 Eg ^ (p)/9p 2 ] p= o, where Eg°\p) is the 
ground-state energy as a function of p for a given value of the dimensionless channel conductance g. (One might 
choose to speak of an effective differential capacitanceQ C c ff = e 2 /2[/ c ff.) SimilarlyEBEj up to a proportionality factor, 
the first derivative I p can be understood as a tunneling-induced correction to the expectation value of h. 

What is desired here, however, is / itself, I being proportional to Diagram 2's contribution to the fourth-order 
energy shift [recall Eq. (A9)]. Integrating I pp twice with respect to p gives us / up to additive terms that are constant 



or linear with respect to p: 



I{p, F, ip) = a + flip + / dx x 



dx,' 



(A17) 



As mentioned in Sec. II and at the beginning of this appendix, the unknown terms (do + dip) are not relevant to 
our result. The a\p term is negligible due to the existence of the mirror image of Diagram 2, in which the roles of 
dots 1 and 2 are exchanged. Such a switch of h\ and fi 2 is equivalent — in calculating energies — to taking p — > —p. 
Consequently, when the total fourth-order shift is calculated, the a\p in Eq. (A17) cancels with the —a\p from the 
mirror image. Likewise, the ao part drops from the final result as we are only concerned with the difference between 
the energy shifts for arbitrary p and p = 0. 

The irrelevance of the ao and a\p terms tells us that we need only calculate / modulo terms constant or linear with 
respect to p. In other words, we need only find an equivalence class 



I(p, F,i/>) - / dxi 



dx2lpp(x 2 ,F,ip) 



(A18) 
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where the congruence symbol indicates equivalence up to additive terms that are constant or linea r wit h respect to 
p. We are therefore free to drop any constant or linear terms that crop up on the right side of Eq. ( |A18| ) . 

Confident that we have figured out wha t we wish to do, we can return to the pedestrian business of doing it. We 
observe that the p-dependent sum in Eq. ( A16) can be written in a more integrable form: 



E 



(-1) 



n+l 



1-P 

3-p 



E 



(-1) 



n+l 



dx In 



In2[ln(l-p)-ln(3-p)+ln3] 
-1 1 



2 - x 



3 — x J \1 — x 3 — x 



Integration of Ip]) and 1$ with respect to p gives 



2 " 



ln(l-p) 



1 °° 

ln([l - F]ip) ]a(F1>) - - ln 2 (F^) + £ 



(_l)n+l 



-t + E 



(_!)«+! /I 



n=l 



n=l 

In 2 In 3 
p ln 2 (3-x) 



+ ~ ln 2 (l - p)[ln([l - F]ip) - ln2] - \ £ dx ±j 



^ln(3-x)ln(l-,)_ ln2 f-J 



dx\ 



p-xi 



l-x 

Xl 



dx2 In 



2-x 2 



•1 



— a; 
ln(3 - x) 
l-x 
1 



1 — Xl Jo " \i — X2 ) \\ — X2 ' i — X2 

for F < (1 - F), 



(A19) 



F 



1 - F 



2' 



ln(l-p) 



1 OO 

-ln 2 ([l -F]^)-]T 



(-1) 



n+l 



1-F 
F 



E 

n=l 



(-1) 



n+l 



In 2 In 3 



[- ln(3-p) + ln3] 
1 



for F > (1 - F); 

(_!)«+! 



6 ^ 



n=l 



- hi21n3 



L ln d (3 - p) - In- 3 3] - - In 2^(3 - p) - hV* 3] 

1 /- p , ln 2 (l-x) T , ln(3 - x) lnfl - a;) 

ax /ax 

3 x Jq 3 x 

/' , lnfl - x) 
ln2 / dx— i ^ 

3 — x 

•X! 



dxi 



3 — xi 



dx 2 In 



-1 



2-x 2 

3 - x 2 J V 1 — x 2 ' 3 - x 2 



(A20) 



The ellipsis in the second equation for Ip indicates that the remainder of 1^ for the system above half-filling is the 
same as the corres pondi ng remainder for the system below or at half-filling. .— . 

In deriving Eq. (A20), we eliminated a number of integrals over x^'s by using an identityE3 that is easily derived 
for double integrals: 



dxi / dx 2 /(x 2 ) 



dx(p-x)/(x). 



(A21) 



Nonetheless, in the final terms of Ip and Ip , double integrals remain. These can be reduced to single-integral form 
with a little extra work. Defining L^™ 1 * 1 to be the last term of (1/2) Jp , we discover that 
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dx2 In ( - 

o V 3 - x 2 J V 1 — x 2 3 - x 2 



■1 1 

+ 



X2 V 



1 — X\ 



- I dx In [ — . | 

V 3 — a; / V 1 — x 3 — x 



o 



-1 1 

+ 



(p - x) + (1 - p) In 



1-P 
1-x 



Jo \3-xJ\l-x 3-xJ 



(p - x) + (3 - p) In 



3-P 

3 - a; 



Now we perform the final integration over p, deriving 1^ and 7^ 2 \ where / = 1^ + /( 2 ). 

1 1 °° f— 1 V'^ 1 / P 1 

(l-p)ln(l-p) ln([l-F]^)ln(F^)--ln 2 (^) + E " ' ' 



6 ' — ' 



2 ^— ' n 2 V 1 - F 

n=l v 

In21n3 + ln([l - F]tp) - In 2 
ln 2 (3-x) 



- \{l P) ln 2 (l - P)[ln([l - m -\n2]- l -j^dx{p- x)^- 

fP . ,ln(3 - x) ln(l - x) , „ , . ln(3 - x) 

+ / <ix(p-x) — i '- — i '--\n2 / dx(p-x) 

Jo 1 - x Jo 

+ (1- P )^ln(l^)ln 



1-x 



2 - x 



-1 1 

+ 



p [2- x 

dx (p — x) In 



3 — x / \ 1 — x 3 — x 

! 1 



3 — x J \1 — x 3 — x 
for F < (1 -F), 



(A22) 
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^ i - +ln21n3 + ln([l-F]V0-ln2 



71=1 
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1/(2)^ [( 3 -p)ln(3-p)-31n3] 
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IT 



71=1 
.2, 



+ - [(3 - p) ln 3 (3 - p) - 3(3 - p) ln 2 (3 - p) + 6(3 - p) ln(3 - p) 

-31n 3 3 + 91n 2 3- 18 In 3] 
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^ (p _ x)h (|^)(_zL + _l_l 



We are essentially done. Upon multiplying the sum of and/( 2 ) by N ch U 2 (t/S) 4 , we have the relevant contribution 
from Diagram 2 to the fourth-order energy shift. After so much work, one might wonder whether we have achieved 
anything more. Providentially, the answer is that, yes, we have. As explained earlier, we have also solved f or th e 
contribution from the corresponding mirror- image diagram, which is obtained by replacing p with — p in Eq. ( |A23| ) . 
Perhaps more surpr isin gly, we have solved for the contributions from another pair of mirror-image terms. A swap of F 
and (1 — F) in Eq. ( |AS| ) turns it into the formula for the contribution from Diagram 3 of Fig. 2(b). Thus, exchanging 



F and (1 — F) in Eq. (A23) yields the contribution from Diagram 3. A further replacement of p with — p gives the 
contribution from Diagram 3's mirror image. The cost of calculating Diagram 2 is high, but at least we benefit from 
a package deal — 4 for the price of 1. 



APPENDIX B: DETAILS OF THE STRONG-COUPLING CALCULATION 



This appendix consists of three parts presenting various calculations described or cited in Sec. III. The first part 
computes \ K c (t), and K c (0), thereby producing the results quoted in Eqs. ( p7| ) and making explicit the 

origin of the factor e 7 that appears in the prefactors of Eqs. (^0|) and (|34|) . The second part of the appendix provides 
the derivation of the first strong-coupling energy correction [see Eq. (|33|)] . The third part derives the second strong- 
coupling correction [see Eq. (]34)]. 



1. Calculation of Charge-Channel Averages 

The leap from Eq. ( |26| ) to Eq. (j27j) in Sec. Ill requires evaluation of the expectation values 

D\ = (cos 



V^c(T) + ^ 



Do 



up 



irp 



(Bl) 



[Recall that time-ordering is implicit in the path- integral definition of (^4),; in Eq. (J24|).] The cosines and products of 
cosines can be written as linear combinations of terms of the form e , where Z is linear in the charge displacement 
operators c (t) and the charge displacement operators arc themselves linear in boson creation and annihilation 
operators (see Ref. 21). Therefore, one can apply a standard relation for the expectation value of the exponential of 
a linear combination of boson operatorsEj 



which can easily be shown to hold for our charge-integration brackets with implicit time-ordering. 
Using Eq. (B2), we discover that 



(B2) 



Di = e~* Ko(0) cos 



D 2 = i e -[^(0)+K c (rx-r 2 )] cos(?rp) + ^[K^-K^-r,)] 

where K c (t) is the charge-charge correlation function, 

K c {t) = (0 c (r)0c(O)) c ■ 

Eq. (27) of Sec. Ill follows immediately. We find S^ 1 ' by replacing cos[v / 7r0 c (r) + 7rp/2] in 5b with D\. For 5 b 2 \ we 
recal 



(B3) 



(B4) 



lat 



([5 b - (5 b ) c ] 2 )c — (5 b ) c - (5 b )c 
and apply the formulas for D\ and D 2 accordingly. 



22 



To get the formula for K c (t) [Eq. (p8|)], we must labor a bit more. Because the unperturbed action Sq^ is quadratic 

charge disp 
Consequently, 



in charge displacement operators c (u> m ), Sq fits exactly the form for the canonical action of a real scalar field 



9 c (u) m )9 c (-uj n )) =- ; S u Un . (B5) 

From this identity and the relation between 9 c (t) and its Fourier transform [recall Eq. Jig) ], we construct a summation 
formula for K c (t): 

In the zero-temperature (/3 — > oo) limit, we may safely transform this sum into an integral. Before doing so, however, 
we should note that, unless the oj m possess an ultraviolet cut-off, K c (0) diverges logarithmically. The standard means 
of imposing such a cut-off in Luttinger liquid theoryc2 is to insert a factor of e^^ m ^ w on the right side of Eq. (|lg|). 
This insertion generates a factor of e - 2 \ u ™\/ w m Eq. (B^), yielding 

r°° du e~ i(AjT e~ i w L 

which is equivalent to Eq. ( |28| ) in Sec. III. 

The way is clear for evaluation of the same-time corrrelation function K c (0). After setting r = in Eq. (B7), we 
integrate by parts and convert to the dimensionless integration variable x — 2ui/W. The result is that 



7T \irW J 7T J 




4U 2 /ttW 



- ln (^7 J +-e (4C/2/ " w) [ / dxe~ x \nx- I dx e~ x In x I . (B8) 



nW ) 7r 

The first integral in the parentheses equals the negative of 7, the Euler-Mascheroni constant.lifl The second integral 
goes to zero as we take the limit W/U2 -^00. In this limit, the exponential factor multiplying the integrals goes to 
1. The final result is the following: 

The derivation of Eq. ( |B9| ) shows that the coefficient e 7 comes from exponentiating a secondary part of (0 C (O)6' C (OU_ 
One might be concerned that Luttinger liquid theory does not faithfully capture such subsidiary dependences. l3 
However, Sec. IV presents evidence that these coefficients are general and independent of the high-energy band 
structure. 



2. The First Strong-Coupling Correction 

As stated in Sec. Ill, in the limit of strong coupling (g — > I), the first correction [see Eq. (p3Q l to the open-channel 
(g = 0) ground-state energy is obtained by diagonalizing the Hamiltonian H New [see Eq. (|29D] . This-,diagonalization 
can be accomplished through another version of the "debosonization" procedure used by Matveev.0 As we wish to 
"debosonize" the action Snow = Sq + [recall Eqs. (|2l]) and J27|)], it is useful to observe-that Sq corresponds to 
the Euclidean action for non-interacting fermions on a semi-infinite lattice ending at x = O.cJ For these fermions, we 
take 9 s {t) to correspond to the x = value of the phase field, <pf(r) = <&/(0,t), rather than the x = value of the 
charge displacement field 9f(r) = 6/(0, t). Making #/(t) = the boundary condition at the edge, we find that the 
properly normalized creation operator for a fermion at x = is given by 
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where, as usual, W is the bandwidth and vf is the Fermi velocity.llB (0, r) can be expressed in terms of reciprocal- 
space creation operators: 

$(0,t) = --= y_ A */I. (Bll) 
The fermionic energies are cut off in the usual way at W/2, the corresponding wave- vector cut-off being A = W/2hvp. 

(s) 

After these machinations, "refermionization" proceeds apace. Since the unperturbed action Sq is an action for 
non-interacting fermions, the unperturbed Hamiltonian Hq is simply the sum of the single-particle energies of those 
fermions. On the other hand, the perturbation that corresponds to is a term linear in fermion creation and 
annihilation operators. In particular, using Eq. (|B9D to determine e~% Ka (°\ we obtain 

= / dk^kftfk, 



^ = V cosf [ A dk (ft + A ) , (B12) 



Not being quadratic in fermion creation and annihilation operators, the. fermionic Hamiltonian we have derived is 
not yet in an easily diagonalizable form. To make it so, we follow MatveevQ in defining a new set of fermion operators 
such that 

fk = [d + S)d k . (B13) 



Plugging this equivalence into Eq. (|B12j) yields Eq. (pCf ) of Sec. III. 

One can now perform the Bogoliubov transformation that produces Eq. (|3l"|). To find the correction to the open- 
channel energy, one notes that H^' of Eq. ( pp| ) has an expectation value of zero in the ground state of , which 
is the open-channel (V — 0) part of H^ evr . Therefore, if the ground state of Hq is represented by the ket |0), 

(0|i?New|0) = Eq, where Eq is the ground-state energy for Hq. From the diagonalized form of H^ cw [see Eq. ([§l|)l, 
it is then deduced that the equation for (E-^ ew — Eq) is the following: 

= - [ A dk^{0\ClC k + & k C k \0) . (B14) 



At this point, it is necessary to know the exact equations for Ck and Ck- As found by Matveevjj for V = 
V 2 [8e^U 2 /TT 2 } cos 2 (7rp/2), they are 

- d k + d)_ k 




n 6c d k - d -k / hv F T j 
C k = 7^=^ 7k I r 2 ) ( + J) 

dy - d)_ k , 



(B15) 



V2 ' 

As before, the symbol V indicates that only the principal value of the integral is computed. 
With the explicit equations for Ck and Ck before us, it is clear that, for k > 0, Ck\0) — 0, and 

AlH(p) = - [ A dkCk(0\CiC k \0). (B16) 



Concentrating on what remains, we see that, for k > 0, both the first term of Ck and the k' > part of the third 

(s) 

term of Ck annihilate the ' ground state. Hence, 



(oictoio) h "' T 



hv F v r 2 (i l 



2^(e' + r 2 ) ^ 2 (£| + r 2 ) \k k + A 



(B17) 
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Plugging into Eq. ( B16| ), we find that 
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47T 



111 



7T 
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Here we have dropped terms that vanish as W/U2 
V 2 [9,eiU 2 /TT 2 ] cos 2 (?rp/2). 
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2^ 




v 2rJ 





W/2 



£, 



f)(^ 2 



r 2 



r 2 ) 



(bis) 



00. Eq. ( |33| ) is obtained by applying the identity T 



3. The Second Strong-Coupling Correction 



The second correction term in the strong-coupling limit [see Eq. (S3)] is derived by treating [see Eq. (5t|)] as 



a perturbation to the system described by H^ cw of Eq. (1 
potential in the finite-temperature path-integral approach/ 



Using the standard formula for the grand-canonical 



f2o = — 4 y^(All connected graphs) , 



we see that the lowest-order correction to the ground state energy of H^ cw is given by 



A$(p)= lim ^(New|^|New>, 



1 



?( 2 )l 



(B19) 



(B20) 



where |New) is the ground-state ket for -£Tncw The minus sign in Eq. (B1E) has been canceled by the minus sign that 



(2) 

arises fcpm the fact that this leading term from St corresponds to a first-order graph and therefore carries a factor 



of-1 



Recalling Eq. d27| ) and observing that the parts of that are independent of p are irrelevant to calculation of 
the fractional peak splitting /, our immediate task is to evaluate the quantity 



*(n,r 2 ) 



vw 



-ttKo(o) cos 2 (^) (New| cos[^ s (ri)] cos[^ s (ri)]|New) 



(B21) 



Under "debosonization" (see Part 2 of this appendix) , this becomes 



-A 



X(n,T 2 ) = A 2 / dfci / dk 2 (New| \dl ± (d + it) + (d + i f )i fel 



-A 



4 (d + d)) + {d+S)d k J |New) : 



(B22) 



where the bracket subscripts indicate that the enclosed operators are evaluated at imaginary times t\ and t 2 , respec- 
tively, and we have used 



A = V cos{irp/2)^2e~<hv F U 2 /iT :i . 



(B23) 



We are now within hailing distance of Eq. (p4[). Using the truncated equations for C k and C k [recall Eq. (|32|)], we 
express the i^'s in terms of these operators. The sub-leading terms in this transformation are negligible as, in the 
end result, they take us beyond second order in V. Similarly, the time-dependence of the operator sum (d + d') is 
sub- leading as (d + cF) first appears in the expansion of the diagonalizing operators at order V. Accordingly, (i+ <v) 
commutes with i?Ncw to zeroth order and can be considered time-independent. In contrast, from Eq. pl|), we know 
that Cfe(r) = Cke~^ kT and C\(t) = C\e^ kT . Application of these insights to Eq. ( B22j ) gives 
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X( n ,T 2 ) = 2X 2 [ dkj rffc 2 (New|C fel (ri)C| 2 (T 2 )|New) 
Jo Jo 



2A 2 



W/2 



hVF Jo 

2A 2 1 - e -(n-T2)W/2 



(B24) 



We now return to Eqs. ( |27| ) and (B20). Switching to dimensionless variables = TiW/2 and substituting for A, we 
obtain 

A (% , f/2 8e^t/ 2 f^ 2 f 1 

KM = V cos \-)—J d *iJ d ^ 



X 1 



--nK c [2(x 1 -x 2 )/W] 



-(xi— x 2 ) 



Xi - X 2 



(B25) 



We eliminate one of the integrations by expressing the integrand in terms of x = [x\ — x 2 ) and observing that in the 
double-integral the density of states for a given value of x is ((3W/2 — x): 



4%) = ^w(f)J 7[/2 



3 PW Jo 



(3W/2 



dx 



-ivK c (2x/W) 



(B26) 



Transformation of Eq. (B26) into Eq. ( p4[ ) follows recognition of the fact that, for x on the order of pW/2, the 
integrand is effectively zero. This is known from the identity 



-Re jV^V^Ki+^Ei [- {AU 2 /nW) (1 + tar)]} 



(B27) 



where Ei[— z] is the first exponential integral function!! For z ^> 1, Ei[— z] goes as e~ z / z. Therefore, the integrand 
goes to zero as I fx 2 for x > irW/4U2, and the region x ^> W/U2 makes a comparatively negligible contribution to the 
integral. This conclusion corroborates the statement made in Sec. Ill that the factor [1 — e~ nK °( 2x > w '] f urrusnes an 
ultraviolet cut-off on the order of ip = W fUi- Since we calculate in the limit — > 00, we know that (3W/2 ^> W/U 2 
and, hence, that the integrand is effectively zero for x on the order of (3W/2. We can approximate the weight function 
(/3W/2 - x) by (/3W/2). The result is Eq. @. 



2G 



1 For an introduction to "single-electronics," see M. A. Kast- 
ner, Rev. Mod. Phys. 64, 849 (1992); D. V. Averin and K. 
K. Likharev, in Mesoscopic Phenomena in Solids, edited 
by B. L. Altshuler, P. A. Lee, and R. A. Webb (North 
Holland, Amsterdam, 1991); and several articles in Single 
Charge Tunneling, Vol. 294 of NATO Advanced Study In- 
stitute Series B Physics, edited by H. Grabert and M. H. 
Devoret (Plenum, New York, 1992). 

2 F. R. Waugh, M. J. Berry, D. J. Mar, R. M. Westervelt, 
K. L. Campman, and A. C. Gossard, Phys. Rev. Lett. 75, 
705 (1995). 

3 F. R. Waugh, M. J. Berry, C. H. Crouch, C. Livermore, D. 
J. Mar, R. M. Westervelt, K. L. Campman, and A. C. Gos- 
sard, Phys. Rev. B 53, 1413 (1996); F. R. Waugh, Ph.D. 
thesis, Harvard University, 1994. 

4 J. M. Golden and B. I. Halperin, Phys. Rev. B 53, 3893 
(1996). 

5 K. A. Matveev, L. I. Glazman, and H. U. Baranger, Phys. 
Rev. B 53, 1034 (1996). 

6 K. A. Matveev, L. I. Glazman, and H. U. Baranger, 



preprint (12/11/95, |cond-mat/9512082[ ). 

7 D. S. Golubev and A. D. Zaikin, Phys. Rev. B 50, 8736 
(1994). 

8 Hermann Grabert, Phys. Rev. B 50, 17 364 (1994); Physica 
B 194-196, 1011-1012 (1994). 

9 K. A. Matveev, Phys. Rev. B 51, 1743 (1995). 

10 L. W. Molenkamp, Karsten Flensberg, and M. Kemerink, 
Phys. Rev. Lett. 75, 4282 (1995). 

11 Karsten Flensberg, Physica B 203, 432 (1994); Phys. Rev. 
B 48, 11 156 (1993). 

12 G. Falci, J. Heins, Gerd Schon, and Gergely T. Zimanyi, 
Physica B 203, 409 (1994); G. Falci, Gerd Schon, and 
Gergely T. Zimanyi, Phys. Rev. Lett. 74, 3257 (1995). 

13 S. V. Panyukov and A. D. Zaikin, Phys. Rev. Lett. 67, 
3168 (1991); Physics Letters A 183, 115 (1993). 

14 H. O. Frota and Karsten Flensberg, Phys. Rev. B 46, 
15 207 (1992). 

15 C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992); Phys. Rev. B 46, 7268 (1992); ibid. 46, 15 233 
(1992). 

16 Gordon Baym, Lectures on Quantum Mechanics (Addison- 
Wesley, Redwood City, 1973), 113ff. 

17 M. P. A. Fisher and Wilhelm Zwerger, Phys. Rev. B 32, 
6190 (1985). 

18 Murray R. Speigel, Mathematical Handbook of Formulas 
and Tables, Schaum's Outline Series (McGraw-Hill, New 
York, 1992). 

19 L. I. Glazman and K. A. Matveev, Zh. Eksp. Teor. Fiz. 
98, 1834 (1990) [Sov. Phys. JETP 71, 1031 (1990)]; K. A. 
Matveev, ibid. 99, 1598 (1991) [72, 892 (1991)] 



20 W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. 
P. Flannery, Numerical Recipes in C: The Art of Scientific 
Computing, 2nd ed. (Cambridge University Press, Cam- 
bridge, 1992), 161. 

21 F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981); J. 
Phys. C 14, 2585 (1981). 

22 V. J. Emery, "Theory of the One-Dimensional Elec- 
tron Gas," in Highly Conducting One-Dimensional Solids, 
edited by J. T. Devreese, R. P. Evrard, and V. E. van Doren 
(Plenum, New York, 1979), 247-303. 

23 I. M. Ruzin, V. Chandrasekhar, E. I. Levin, and L. I. Glaz- 
man, Phys. Rev. B 45, 13 469 (1992); L. I. Glazman and 
V. Chandrasekhar, Europhys. Lett. 19, 623 (1992). 

24 C.J.B. Ford, P.J. Simpson, M. Pepper, D. Kern, J.E.F. 
Frost, DA. Ritchie, and G.A.C. Jones, Nanostructured 
Materials 3, 283 (1993). 

25 C. A. Stafford and S. Das Sarma, Phys. Rev. Lett. 72, 3590 
(1994). 



C. A. Stafford and S. Das Sarma, preprint (5/12/95, cond 



mat/9505058) 



77 M. Kemerink and LW. Molenkamp, Appl. Phys. Lett. 65, 
1012 (1994). 

28 N. Tsukada, M. Gotoda, M. Nunoshita, Phys. Rev. B 50, 
5764 (1994). 

29 T. Sakamoto, SW. Hwang, F. Nihey, Y. Nakamura, and 
K. Nakamura, Jpn. J. Appl. Phys. 33, 4876 (1994); Super- 
lattices and Microstructures 16, 291 (1994). 

30 G. Klimeck, Guanlong Chen, and S. Datta, Phys. Rev. B 
50, 2316 (1994); Guanlong Chen, G. Klimeck, S. Datta, 
Guanha Chen, and W. A. Goddard Iff, ibid. 50, 8035 
(1994). 

31 F. Hofmann, T. Heinzel, DA. Wharam, J.P. Kotthaus, G. 
Bohm, W. Klein, G. Trankle, and G. Weimann, Phys. Rev. 
B 51, 13 872 (1995). 

32 R.J. Haug, R.H. Blick, and T. Schmidt, Physica B 212, 
207 (1995). 

33 N.C. van der Vaart, S.F. Godijn, Y.V. Nazarov, C.J.P.M. 
Harmans, J.E. Mooij, L.W. Molenkamp, C.T. Foxon, Phys. 
Rev. Lett. 74, 4702 (1995). 

34 C. H. Crouch, C. Livermore, F. R. Waugh, R. M. West- 
ervelt, K. L. Campman, and A. C. Gossard, Proceedings 
of Electronic Properties of Two-Dimensional Systems XI 

Surf. Sci. (to be published)] 



35 P. Pals and A. MacKinnon, pre print (12/31/96, |cond- 
mat/96011511 ; preprint (12/31/96, |cond-mat/9601156^ 



D. J. Amit, Field Theory, the Renormalization Group, and 
Critical Phenomena, 2nd ed. (World Scientific, Singapore, 
1984). 

37 J. Solyom, Adv. Phys. 28, 201 (1979). 

38 J. W. Negele and H. Orland, Quantum Many-Particle Sys- 
tems (Addison- Wesley, Redwood City, 1988). 

39 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Se- 
ries, and Products, 5th ed., edited by Alan Jeffrey (Aca- 
demic, Boston, 1994), 933ff. 



27 



